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TUNNEL , WITH CONSIDERATION 0? THE 
EFFECT OF COMPRESSIBILITY 
By H. Julian Allen and Walter G. Vincenti 

■ SUMMARY 

Theoretical tunnel-wall corrections are derived for an air- 
foil of finite thickness and camber in a two-dimensional-flow 
wind tunnel. The theory takes account of the effects of the 
wake of the airfoil and of the compressibility of the fluid, and 
is. based upon the assumption that the chord of the airfoil is 
small irl comparison with the height of the tunnel. Consideration 
is given to the phenomenon of choking at high speeds and its 
relation to the tuimel-vall corrections. The theoretical results 
are compared with the small amount of low- speed experimental data 
available and the agreement is seen to be satisfactory, even for . 
relatively large values of the chord-height ratio. 



HJTRODUCTION 

The need for reliable wind-tunnel data for the design of 
high-performance aircraft has led in recent years to attempts to 
make the conditions of • the 'tunnel, tests conform more closely with 
the conditions prevailing in flight, especially with regard to 
the Reynolds and Mach numbers. Because of practical limitations 
in size and power, most existing wind tunnels, whether high speed 
or low speed, are not capable of providing full-scale Reynolds 
numbers for all flight conditions. In order to obtain the 
highest Reynolds numbers possible under the circumstances, it. is 
necessary to use models dimensions of which are as large as 
possible relative to the cross-sectional dimensions of the tunnel 
test section. The effect of such large size is to make the test 
conditions -^epa^t further from the conditions prevailing in 
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flight by increasing the magnitude of the tunnel-wall interfer- 
ence. In the case of teet.i at high Mach number 8, the interfer- 
ence is increased still further by the tendency of the flow 
pattern of a compressible fluid, if unrestrained, to expand as 
the Mach number of the undisturbed stream increases. Since .the 
walls of a closed-throat tunnel restrain certain of the stream- 
lines at a fixed distance from the model, this expansion is 
prevented; and the taimel-wall : interference and corrections be- 
come progressively larger as the Mach number increases. The 
results obtained in the tunnel must therefore be corrected accu- 
rately for the effects of. wall interference if they are to be 
applied with confidence to the prediction of free-flight charac- 
teristics. 

In tests at high Mach numbers an additional complication 
arises. The effect of a model in. a closed-throat tunnel may, in 
a sense, be thought of as equivalent to that of a constriction 
in the throat of the tunnel. The resulting converging-diverging 
nozzle formed by the model and the tunnel walls then has roughly 
the same characteristics at high speeds as the usual supersonic 
nozzle- that is, for some Mach number less than unity in the 
undisturbed stream, sonic velocity is reached at all points 
across a section of the tunnel at the position of the model, and 
the flow in the diverging region downstream of this section be- 
comes supersonic. When this occurs, increased power input to 
the tunnel has no effect upon the velocity of the stream ahead of 
the model, the additional power serving merely to increase the 
extent of the supersonic region in the vicinity of the model. At 
this point the tunnel, is- said to b.e "choked" and no further in- 
crease in the test Mach number can be obtained. The value of the 
Mach number at -which choking occurs is thus of extreme importance, 
since it determines the upper limit of the range of Mach numbers 
which can be obtained with a given combination of model and 
tunnel. 

.In testing airfoils to obtain section characteristics at . 
subsonic speeds, it has become common practice in modern closed- 
throat wind tunnels to have the model 3pan the tunnel so that 
supporting struts and their accompanying interference effects are 
entirely eliminated. If the tunnel has a cross section of 
rectangular shape, this arrangement results in a flow which is 
essentially two-dimensional. 

The wall interference for such a two-dimensional-flow wind 
tunnel has been the subject of numerous investigations, the results . 
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in general being expressible as series in ascending powers of 
(c/h) , where c is the chord of the airfoil and h the height 
of the tunnel. The effect of wall interference upon the flew of 
an ideal fluid about a symmetrical airfoil at zero angle of attack 
is determined to the order (c/h) 3 by Lock in reference 1 and 
by Glauert in reference 2. The interference for an infinites!- 
mally thin, cambered airfoil at a small angle of attack in an 
ideal fluid is given by Glauert to the* order (c/h) 2 in reference 
2, and investigations for the special case of a flat plate have 
been carried out to a higher order of accuracy by several writers. 
While the present report was being prepared, work by Goldstein 
appeared (references 3 and 4) in which the interference is 
determined to the order (c/h) 4 for a general cambered airfoil 
of finite thickness in an incompressible fluid, no restriction 
being made in the general results as to the magnitude of the 
camber, thickness, and the force coefficients. A still later 
paper by Goldstein and Young (reference 5) gives the modifications 
necessary in the previous results to allow for the effect of 
fluid compressibility to the order (c/h) 3 . 

In the present paper, the tunnel-wall corrections are deter- 
mined to the order (c/h) s for the general airfoil in a compres- 
sible fluid for Mach numbers below that at which choking occurs. 
It is assumed that the thickness and camber of the airfoil are 
small and that the interference velocities are everywhere small 
as compared with the velocity of the undisturbed stream. A dis- 
cussion is also included, of *the Mach number at which choking 
occurs. The various results presented are of essentially the 
same nature as those which already h$vo appeared separately in 
the reference cited, but the methods of development and certain of 
the final results, are different, especially with regard to the 
interference associated with lift. The validity of the final 
corrections is examined by comparison with the available experi- 
mental data. The equations also are compared with the results' of 
references 3, 4, and 5, and the aforo-mentioned differences are 
discussed. 

The discussion is limited to airfoils placed midway between 
the upper and lower walls of the tunnel. Mathematical symbols 
are defined as introduced in the text. For reference, a list of 
the more Important symbols and their definitions is given in 
appendix B. v 
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Wm/SFOSm 0? COPMCTIOE EQUATIONS 

" In an analysis of tunnel-vail interference it is desirable 
to look upon the theoretical development of the tunnel-wall 
corrections as consisting of two parts. First, it is necessary 
to determine the rnanner and extent to which the tunnel walls 
alter the field of flow about the airfoil from what it would be 
if they' were not present/ Second, it is necessary to calculate 
the effect of thase alterations upon the measured characteristics 
of the airfoil. The development of the correction equations of 
this report has been divided, into those two general sections. 

In reference 6, the use of the method of superposition to. 
determine the pressure distribution over the surface of an air- 
foil section in free air is presented. It is shown that in the 
calculation cf the flow .at the surface of a thin airfoil of small 
camber, the effects. of. camper and thickness may be considered 

^'independently. This- .follows directly from the fact that the 
velocities induced by the vortex sheet used to represent camber 
and those induced by the source-sink system used to represent 
thickness are simply additive in their effect on the flow over 

' the airfoil . 

To treat the problem of wall interference, it is. again 
convenient to consider . the thickness and camber effects separately. 
The flow changes associated with airfoil thickness are found by 
considering the interaction' between the tunnel wails and the base 
profile of the airfoil, the base profile being defined as the 
profile the airfoil would have if the camber were removed and the 
resulting symmetrical airfoil placed' at. zero angle of attack. 
The interference effects associated with airfoil camber are found 
by analyzing the interaction between the' tunnel walls and an 
' infinites irnelly thin airfoil having the same camber as the actual 
airfoil. In addition to'the interference effects associated with 
airfoil thickness and camber, it is necessary to consider a fur- 
ther alteration of the field of flow caused by the confining 
influence of the tunnel walls, upon the airfoil wake. When the 
individual effects promoted by the interference between the walls 
and the airfoil thickness, camber, and wake are known, the total 
alteration in the flow at the airfoil is found by .superposition, 
and the characteristics of the airfoil in the altered field of 
flow are compared with the characteristics in free air. This 
comparison. leads to simple formulas which enable the prediction of 
the free- flight characteristics when the characteristics in the 
tunnel are known. 
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The method, of superposition, which is f undamental to the 
entire analysis, la in general inapplicable to compressible flow • 
as the differential equation for such flow is nonlinear in the 
physical plane. The separate solutions which are superposed are 
obtained, however, by assuming that the airfoil is of small thick- 
ness and camber and that the induced velocities are thus small ai 
compared with the velocity of the undisturbed flow. On the basis 
of this assumption the equation of compressible flow becomes a 
linear differential equation - namely, Laplace 1 3 equation (refer- 
ences 7, 8 ? and 9) - so that superposition of velocities is, in 
this case, technically permissible. . Furthermore,- the tunnel-wall 
corrections are in most cases rather small relative to the experi- 
mental quantities being corrected, so that .it is not thought -that 
the use of this approximate method will lead to large errors in 
the final corrected quantities. 

Influence of Tunnel Walls upon Field 

of Flow at Airfoil 

Thickness effect.- The interaction between the base profile 
and the~walis °* a ""two-dimensional-flow tunnel has been consider- 
ed by Lock for the case of an incompressible fluid (reference 1; 
a discussion of Lock's method is also given by Glauert in 
reference 2)". Lock's method of analysis is essentially to intro- 
duce an infinite series of images of the base profile such as to 
satisfy the condition that there is no flow normal to the walls, 
to replace each image by a suitable source- sink' doublet, and to 
calculate the 'velocity induced at the base profile by this 
system of doublets. It is shown that the net effect of the 
tunnel walls upon the flow at the base profile is to increase 
the effective axial velocity of an incompressible stream by : the 
amount t 

(4 V), ssAOV' '(1) 

where 

V 1 apparent stream velocity at airfoil as determined from 
measurements taken at a point far ahead of model 

a f. factor dependent upon size of airfoil relative to tunnel 

A a factor dependent upon shape of base profile 

The factor O is. defined by the equation 
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where ' (c/h) is the ratio of the airfoil -.chord to the tunnel 
height. The factor A can he determined for any '^base Drofile 
from the relation 



A = Ifif £fc 




2 - 



(!) 



(3) 



where 



ordinate of base profile at chordvise station x 



fiy-t/dx slope of surface of "base profile at x 

7f. base-profile pressure coefficient at x 
1 compressible fluid 



m an m- 



(It will be noted that the quantity Xtr ' in references 1 and 2 

is equivalent to -iAc in the notation of this report.) Values 
4 

of A for a number of base profiles are given in table I. 

In appendix A, it is shown that the effect of compress- 
ibility upon the stroamwise induced velocity at a given point a 
large distance above or below a body in a uniform stream is such\ 
as to multiply the velocity increment for incompressible flow by 

the factor l/Q. - M 2 ] 0//Q where M ■ is the Mach number of the 
flow far upstream from the body. Applying this result to the 
velocity induced at the base profile by each of the airfoils *in 
Lock 1 3 system of images , it can be seen that for a compressible 
fluid the increase in the effective axial velocity in the tunnel 
Is 



[1 - (MO 2 ] 



3/i 



AcV 



(4) 



where W is the apparent Mach number - that is, the Mach 
number corresponding to the velocity V 1 . 

It should be. noted that equation (4) does not agree with the 
result given in reference 10, in which it is stated that the 
velocity increment in the incompressible fluid should be multi- 
plied by the factor' l/[l - (M 1 ) 2 j to allow for the effect of 
fluid compressibility. ' k critical review of this latter report by 
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it.3 author and by others has disclosed that the result given there- 
in is incorrect. The error in the analysis arose from a failure 
to preserve the shape of the body for which the compressibility 
factor being determined , so that the velocity at the surface 
of the body was influenced by a change in shape as well as by the 
tunnel-wall interference. This difficulty does not arise in the 
analysis of the present report. The result of equation (4) has 
also been obtained by an independent procedure in reference 5. 

Consideration of the symmetry of the base profile and of 
the system of images uaefl by Lock to simulate the effects of the 
tunnel vails indicates that the interaction between the walls 
and the baee profile does not induce velocities- normal to the 
center line cf the tunnel . Similarly the base profile does not 
effect the longitudinal velocity gradient in the tunnel at the 
position of the airfoil. 

Wake effect, - In the wake of an "airfoil moving through a 
real fluid, the~total head cf the fluid is leas than in the 
region outside the wake. This reduction arises from the increase 
in thermal energy caused by fluid friction in the boundary layer 
and in the wake itself and by any shock waves which may exist in 
the vicinity of the airfoil. Considering a section normal to the 
wake , it may be said that the static pressure across the stream 
is nearly constant if the section taken is not too close to the' 
trailing edge of the airfoil. It follows that the reduction in 
total head which exists within the wake must appear almost entire- 
ly as a decrease in the local dynamic pressure cf the fluid. 
This decrease arises primarily from, a reduction in the local 
velocity and secondarily from the reduction in local density which ■ 
accompanies the increased temperature within the wake. Thus, 
since the local velocity and density within the wake are both less 
than in the external flow, the mass-flow rate per unit area is 
less inside the wake than outside. Thjs condition prevails both 
in the tunnel and in free air. In the tunnel, however, the require- 
ment of continuity , of flow between a transverse section upstream 
from the airfoil and a lection across the wake necessitates, in . 
addition, that the mass- flow rate per unit area outside of the wake 
is greater than the mass- flow rate per unit area ahead of the air- 
foil, In order to satisfy this requirement, the velocity in the 
tunnel outside of the wake must be greater than that of the 'in- 
disturbed stream. This fact implies that as the flow proceeds' ■ 
down the tunnel the velocity of the main portion of the stream 
undergoes a gradual increase from the value prevailing in the 
undisturbed stream ahead of the model to some higher value down- 
stream of the airfoil. This does not hold true in free air, where 
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the velocities of the rcain flow upstream and downstream of the 
model are equal. The interference between the wake and the 
tunnel walls thus gives riss at the position of the model to a 
velocity increment and a velocity gradient which are not present 
in an r^ni raited 3tream. Further 9 as required "by Bernoulli's 
equation, the velocity gradient io accompanied by a longitudinal 
pressure grfcdienx which likewise would not exist in free air. 

To determine the magnitude of these effects the procedure 
is briefly as "olicws: "Two stations in the tunnel are considered, 
one far upstream froaa the model and one far enough downstream so 
that the wake has spread to /ohe walls and the velocity is again, 
unifoCT across the tunnel. The difference in static pressure be- 
tween these two stations is evaluated as a function of the measur- 
ed drag of the airfoil-. The pressure gradient at the airfoil can 
be related to this pressure difference and hence to the drag of 
the airfoil by a convenient analytic device, which is essentially 
the same as that U3ed by Goldstein (reference 3). The airfoil 
and. its wake are considered, to be replaced by a fluid source 
located at the position of the airfoil. It is specified that 
conditions far upstream in the resulting hypothetical flow must 
be the same as those existing in the actual stream, With this 
provision, the magnitude of the velocity and static pressure far 
downstream can be determined as functions of the upstream condi- 
tions and the strength of the source . The strength is then re- 
lated to the drag of the airfoil by requiring that the static 
pressure difference promoted between the • two. stations in the 
tunnel by the source flow is the same as that which actually 
exists when the airfoil and wake are present. The tunnel walls 
.can then be replaced by an infinite system of such sources 
directly above &nd below the position of the airfoil at intervals" 
equal to the height of the tunnel. The system of image sources 
alone, however, would induce a small finite 'negative velocity at 
infinity upstream, so that 'it i3 necessary to superpose on the 
flow field an additional uniform flow of equal, velocity in the 
positive direction in order to satisfy the original requirement 
that the conditions far upstream shall be unchanged. The velocity 
of this flow, which is readily determined as a function of the 
scfurce strength and hence of the airfoil drag, then gives the 
velocity increment caused at the airfoil by the interference be- 
tween the wake and the walls, The longitudinal velocity and 
pressure gradients at the position cf the 'airfoil are found in 
terms of the drag by evaluating the. flow induced at that point 
by the image sources. It is apparent that this entire method of 
analysis fails to satisfy the actual condition as regards the 
velocity at infinity downstream. This discrepancy arises out of 
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the fundamental difference "between the actual flow in the wake and 
the source flow by which it is represented and is unavoidable as 
long as this representation is used. 

Consider the flow in a closed two-dimensional-flow wind 
tunnel, as shown in figure 1, At a station far upstream, the 
effect of the model upon the flow is negligible, so that the 
velocity V, the density p', the static pressure p', and 
the absolute temperature T 1 are constant across the stream. At 
a station far downstream, where the wake has spread to the walls , 
the velocity V", the density p", the pressure p", and the 
absolute temperature T" are again constant across the stream. 

The difference between the pressures p 1 and p" can be 
related to the measured drag of the airfoil by means of the con- 
ditions of continuity, conservation of energy, and impulse and 
momentum, together with the state relations for a perfect gas. 
The. condition of continuity is given by 

pt v> = p* V" (5) 

and, if it is assumed that the flow is an adiabatic process, 
conservation of total energy requires that 

SIX. - 112. = gjo (T» - T.) 

2 2 P 

or 

whore 

g gravitational acceleration 

J mechanical equivalent of heat 

Cp specific heat of gas at constant pressure 

In modern wind tunnels the walls of the test section are flared 
slightly to compensate for the growth of the boundary layer on 
the walls, and only the drag of the airfoil therefore need be 
considered. The impulse-momentum equation can be written 

2? = p« - p" + P'(vf - p"(v") 3 

h 
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or 



P'(V') 3 



+ 1 



91 ( L 
p' \v« 



(7) 



where D' is the drag of the airfoil and c d ' the drag co- 
efficient referred to the apparent dynamic pressure q>', < 

The velocity of sound V c ' in the undisturbed stream is 
related' to' the absolute temperature by 

(V) 2 = KT ' = (Y - l)gJ CT ,T' 



(8) 



where y is the ratio of specific heats and E is the gas con- 
stant. By means of this relation, equation (6) can be written 



V" _ 
7T = <1 



L 



* (g - 0- [ 

(y - i)(M'f 



! l/a 



O) 



and, from equation (5), 



-1/2 



P- =< ! 1 



(Y - 1)(M')? 



(10) 



The state equation for a perfect ga3 then .provides the relations 

Vl/2 

(ID 



p P"T n T ,f 
p 7 = p 7 "? 7 = T 7 S 1 



and 



(Y - 1)(M') Z 



P* • 



MM 



(v c «) 3 



p'(V')* (v) a y(v) 2 Y(M') 3 



(12) 
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Substitution IVcin equations (3), (10) > (11), and (12) into 
equation (?) gives 



r 



2 d \h/ v'\t:^. 2 T* 



W 7 (HO 

r 



2 1 r - 1) 

?Y-l)(M') a 



j (Y - 1)(M0 2 



1/3 



from which it can "be found that 



j (V - 1)(M') 2 [ • (7 + 1)(M') 2 

>- J 



J 



/(> - (M') S ] S -Y(M') 2 (I) f 2 +Y(M0 a U - ^ (-2)" 



(13) 



(V + 1)(M0 2 

For airfoils usually employed , the factor c d J ^\ is small. 
Expanding the above expression and neglecting terms containing 

c d ! ^Sj to powers higher than the first gives 

1 1 _ i(iL^L j ml + 7 i*: /ex _iM_a 3 „ . (u) 

(Y - l)(M') 2 | 2 Uy 1 - (M') 2 

Ey means of this relation, together with equations (9) and (10) 
the static pressure difference in terms of q' is obtained from 
equation (7) as 



P' - P" 

= C; 

q.' 



r Y{j4«) 2 i 



/C\ I Y M' , « 

d t ■- \ ! i + — ! — i K 15 

d U/ l 1 - (m s ) 2 f 
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How. consider the airfoil and valce to.be removed and re- 
placed by a fluid source of strength Q. If the flow" con- 
ditions far upstream are maintained unchanged , the macs flow 
fr.r dovnatream is then • 



hPa'Vg" = hp'V + Q 

or 



P < 

p ! V 



= 1 + 



hp'7 1 



(16) 



where the subscript 3 denotes conditions new prevailing a J , 
the latter s cation. For reversible adiabatic flow 



" 7.1 

a - — (mo 3 ! -i-i - 1 > 

2 L\V' / J J 



Since it is to be expected that (V B "/V T ) will -be close to 
unity the right-hand side of this equation may be expanded 

[7vY 1 

in ascending powers of ! [ — - j » 1 ; } and terms contain- 



ing powers higher than the first neglected. Thus, 



P 
P 1 



5- - 1 . ( M ') ■ " 



(17) 



and equation (16) becomes 
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It is reasonable to assume- that the. ratio — - — is small -as 

hp'V 

compared with unity. This solution of the preceding equation to 



the first order in 



Q 



hp'V 
/V 



is 



(14 x L_-JL. 

VV / l _ ( M i) 2 hp, V i 



(18) 



Bernoulli's equation for reversible adiabatic flow can be written 



P' - P H " 



7(M')' 



r 



1 - a - (M') S 

L 2 



''VV 



I" 

\V 



1*1 



. -J J 



v " 
Since ' 8 



[~~ J is close unity, this may 1 'eplaced by- the 
approximate relation 



P' - P E 



q* 



V J ' 



Substitution from equation (18) and neglect of "the term involving 

the square of — - — then gives 
hp'V ■ 



P' 



2Q 



q' 



hp'V 1 - (M«) 3 



(19) 



Comparison of equations (15) and (19) shows- that the source 
strength required to promote the same pressure difference as 
actually arises from the confinement of the airfoil wake is 



Q = 



p'Vc^'c f 



i-U - ( 7 - i)(m'H 

L J 



(20) 
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- The tunnel vails are now replaced "by .an infinite system of 
sources of strength Q. a#vcel h ditrbaiice apart and located 
directly above and below the position of "the airfoil as shown in 
figure 1. 'This image system together with the source which has 
been placed at the position of the airfoil satisfies the require- 
ment that the flow a.t the plane of the tunnel wall shall be 
tangential to the wall. 

As shown in the first of equations (All) of appendix A, a 
source of strength " Q in a unifbrm flow of compressible fluid 
will induce at a distance r from itself a streamwise velocity 



Zrtpr j ✓T~- ,M 3 -(1 - M 2 sin^cp) 



where cp is the polar angle of the point in question and p 
and M are the density and Mach number of the undisturbed stream. 
By virtue of this relation, the streamwise velocity AgV 1 in- 
duced at a point of the s center line of the tunnel by the entire 
system of image, sources is 



. . y r - - - , CQ ^v -. 

^ *p'r j/l - (M») 3 fl ~ (M ! ) 2 sin 
L L 




where r m and cp m are the radial ' distance and the polar angle 

of the point relative to the source a distance mh above or 
below, the center line and p 1 and M 1 are the density and Mach 
number of the undisturbed flow in the tunnel. If the distance 
from the position of the airfoil to the point on the center line 
is denoted by z (taken positive downstream), this equation can 
be written 



A 2 V« = 1 
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or 





coth 



nx 1 



np'/l *-(M')2 2hyi-(M') 



hyi-(M f F 2x 



It can be seen by- setting x = -oc in equation (21) 
that the image sources induce at an Infinite distance up- 
stream a velocity 



In order to satisfy the original requirement that conditions 
far upstream remain unchanged, this velocity must be counter- 
balanced by the superposition of a uniform flow of equal 
magnitude but opposite sign. The addition of this flow at 
all points in the field will, result in a speeding up of the 
general flow at the postiion of the airfoil by the amount 



or, substituting the source strength from equation (20), 




Q 



2/T 

2p'h | 1 - (M l ) J 




Q 



2p'h (M l ) 




j 1 + ( Y - 1)(M') 
L 1-(M«) 2 



2 




. If the factor t is defined as 




The velocity increment induced at the position of the air- 
foil by the interference between the wake and the walls may 
thus finally be written for air (Y = 1.4) as 
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The longitudinal velocity gradient produced at the position 
of the airfoil by the flow; from the image sources can "be found by 
differentiating equation (21) with respect to x and then 
setting x = -0. This gives finally 



dV f d(Aj>Y 7 ) *Q 

= ~ ~T3 /2 

dx cb: 6p*h 2 [l - (M'fj - 

or, by virtue of equation (20) ; 

£11 - f 1 + jZ : ^)(M f ) 8 1 jCC d ?Y,a (24) 
dx "V [1 - (MO 3 ] 3 / 2 ' J 12h 3 

It already hag been noted -that- the interference associated 
with' the thickness of the airfoil has no effect upon the longi- 
tudinal velocity gradient at the position of the' airfoil. It 
will be seen later that this also is true of the interference 
associated with airfoil camber. Equation (24) thus gives the 
total velocity gradient for the complete airfoil and wake. The 
total pressure gradient at the position of the model then is 
given by Bernoulli's equation as 



dx w dx 



or, substituting from equation (24) and setting 7 = 1.4, 



dp = J 1 4 0.4(M') : 



dx , lb 



6h 2 



(25) 



It is apparent from the symmetry of the system of image 
sources that at the center line of the tunnel the interference 
between the wake and the walls ha3 no effect upon the velocity , . 
normal to the direction of the stream. 

It is shewn later in this report that the camber of the air- 
foil does not affect the stream velocity at the airfoil. Equations 
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(4) and (23) together thus give the total increase in velocity • 
for the complete airfoil &M wake. The effective or true 
velocity V at the model is therefore 

r ^ 

V = Y<\ 1 + I . Aa+ L^9Jj^)l rc d 4 (26) 

L [l - (Mt) 3 ] 3/a 1 - (MO 2 J 



It is evident that .a correction to the apparent velocity in 
a compressible flow implies corrections also to the apparent 
density j dynamic pressure, Reynolds number , and Mach number. 
These corrections are readily obtained on the basis of the usual 
assumption that the flow is a&iabatic. It is assumed that the 
correction terms are small as compared with unity, so that squares 
and products of these terms may be neglected. 

The true density p at the model is connected with the 
apparent density p 1 by the insentropic relation / 



T] r-i 

- 11 (27) 



, p , fi . 21^ < M! H (Xf 

2 • I U7 



Substitution from equation (26) gives, after expansion as an 
ascending power series and neglect of correction terms higher 
than the first order, 

(Mil! -At. (M') a & + 0.4(MjQ T 1 (2Q) 



p = p" 1 - 



[l - (M') S ] 3/S 1 - (M«) 2 "j 



The true dynamic pressure q = ipY -3 is related to the 

2 

apparent dynamic pressure q r by the equation 



q = v (SL) (If 
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By means of equations (26) and (28) this can be written to the 
first order as 



L [l-(M') 2 ] 3/S . 1-(M<) 2 d J 



(29) 



The true Reynolds number R is given in terms of the 
apparent Reynolds number E ' by the equation 



B 



where u and u 1 are the coefficients of viscosity correspond- 
ing to V and V, According to Von Karman and Tsien (reference 
11), tho coefficients of viscosity are related to the correspond- 
ing absolute temperatures by 



JL _ ( J\°' 7i 



For reversible, adiabatic flow it can be shown that 



T = T' { 



7 - 1 ,. 



(M<) 3 



7 V 



which after substitution from equation (26) becomes to the 
first ©rder for air (7 = 1.4) 



=T« jl - 0.4(M') 
L i>(M') 3 ] 



Ao - 0.4(M') a [.U0.4(M') 3 '3 TC , \ 
> Wls/s i . (m« )2 d J 



(30) 



(31) 



By means of these relations together with equations (26) and (28), 
the true Eeynolds number may be written 



E=E 1 



[i + J^Z 



. ji , _l-0.7(M') 2 A ? x [l- 0.7(M') 2 ][l + 0.4(M') 2 ] ^ 



2l3 / 



) 8 3 



1 - (M') 2 



(32) 
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The true Mach number M is related to the apparent Mach 
number M' by the equation 

where V„ and V * are the velocities of sound corresponding 

to V, and V*. Since the velocity of sound in a gas is 
directly proportional to the square root of the absolute temper- 
ature alone, this equation may also he written 

-"■'(S)tfir 

With the aid of equations (26) and (31) the true Mach number then 
may 'be written to the first order 

f r - - ^ 

M = M< J i . 1+0,2(M') 3 , q A LlH<),2(M0 8 JU^^Mlgj TB • j (35) 

L r>(M0 2 j 3/2 " 1 - (mo 2 j 

At low Mach numbers, the terms containing Tq' in the 
correction equations are usually negligible as compared with the 
terms containing Ag. At supercritical Mach numbers, however, 
where the drag coefficient is very large, the terms with Tc d J ' 
are predominant, 

Numerical values of the compressibility factor appearing. in 
equations (26), (29) , (32), and (33) are given in table II. . 

Camber effe ct,- The theory of the inf initesiraally thin, 
cambered airfoil in free air is developed by Glauert in reference 
12 (pp. 87-93). In this development the camber line is replaced 
by a sheet of continuously distributed, bound vortices. The flow 
induced at any point on the camber line by this system of vortices 
is obtained by integration and is combined vectorially with the 
flow of the undisturbed stream to give the direction of the 
resultant flow. The distribution of vorticity is then determined 
from- the condition that the resultant flow at all points on the 
camber line must be tangential to the camber line. 
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In the actual calculation of the induced velocity, it is 
assumed that the vortices r.isy "be distributed along the chord line 
rather than along the camber line and that the induced velocity 
at any chordwise station on the camber line is the same as the 
induced velocity on the chord line at the same station. If the 
origin of coordinates is taken at the leading edge of the airfoil 
(fig. 2), with the positive x-axis along the chord line and the 
positive y-axis directed upward, the induced velocity (v) x in 
an incompressible fluid at any point x Q on the chord line is 

"Tax 



o 

where dF/dx is the 'vorti city per unit length at the point x 
and c i3 the chord of the airfoil. The direction of this 
velocity is normal to the x-axis. 

Glauert (reference 7) has shown that a first approximation 
to the velocity induced at any point by a "simple vortex in a 
compressible stream can be obtained by simply multiplying the 
velocity induced at the same point in an incompressible stream by 
the factor 



1 - M 2 sin 2 cp 



where M is the Mach number of the undisturbed flow and 9 the 
polar angle of the point in question as measured from the direction 
of flow of the undisturbed stream. For points on the chord of an 
airfoil which i3 inclined at a small angle to the direction of the 
undisturbed stream the polar angle 9 is small, and the factor 
(35) is sensibly equal to 



A/1 - M 3 



If it is assumed that the. effect of a vortex sheet in a com- 
pressible fluid may be obtained by superposing the effects of 
elementary vortices, the velocity induced at any point x 0 on 
the chord line in a compressible fluid is 



<v = a/i - M / dx (36) 



2« J x - x 0 
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If the undisturbed velocity of the free stream is taken 
equal to the true velocity 'V at the airfoil in the tunnel , the 
condition that the resultant flow shall be tangential- to the 
camber line requires that, for all points on the airfoil, 

(37) 

V dx 

where dy c /dx is the slope of the camber line at x 0 and', a is 
the true angle of attack; that is/ the angle the undisturbed 
stream makes with- the chord line in free air. (See fig. 2.) 

The problem of the infinites iroally thin, cambered airfoil in 
a two-dimensional- flow tunnel can be investigated • by the method 
of images; that is, the effect of the upper and lower walls of 
the tunnel can be simulated by introducing an infinite lattice of 
alternately inverted but otherwise identical image airfoils above 
and below the original airfoil, as shown in figure 3(a). By this 
artifice the direction of flow at the position of the upper and . 
lower walls can be made to coincide with the plane of the walls, 
which is the required condition of flow. As in Glauert*3 analysis 
of the airfoil in free air, the camber line of the. airfoil and of 
each of its images is 'replaced, by a sheet of continuously distri- 
buted vortices, the vortex distribution of all sheets being 
identical in magnitude but alternately reversed in sign. The 
flow induced at any point of the camber line of the original air- 
foil by the entire vortex system is then obtained by integration. 
As before, the distribution of vorticity must be determined so 
that the resultant of the induced velocity and the stream velocity 
is tangential to the camber line of the airfoil. 

For the detailed calculation, the coordinate system is taken 
as shown in figure 3(b). The origin of coordinates is taken on 
the center line of the tunnel at the leading edge of the airfoil. 
The positive x-axis .extends downstream parallel to the undisturb- 
ed flow, and the positive y-axis is directed upwards. It is 
assumed that the vortices may be distributed along the x-axis and 
the induced velocities calculated at points on this axis. This 
arrangement is somewhat different from the system employed for 
the airfoil in free air, where the x-axis was taken along the 
chord line; however, since the angle of attack is assumed to be 
small, the difference is of no consequence. 



22 



NACA ARB No. 4K03 



It is evident from figure 3(b) that, for an airfoil midway 
"between the upper "and lower walls of the tunnel, the axial 
velocity induced at any point on the x-axis by any one image is 
nullified by the velocity induced by the corresponding mage on 
the opposite side of the tunnel. It follows that airfoil ■ camber 
does not affect either the true axial velocity or the longitudinal 
pressure gradient in the tunnel at the position of the model. 

The vertical velocities induced at any point on the x-axis 
by any one image and its counterpart are, however, additive. 
Thus, for corresponding images situated at mh and -mh, 
respectively, the vertical velocity (v^)^ induced at the 

point x Q ■ in an incompressible fluid is 



/ - — sin \q> - ~ i dx 



(v m )i - 2 (-D m 7 d * m 2 



atr m 



or 



/ |L ( x - x D ) dx 

(v« ). =1 (-1)*/ °J- (38) 

m 1 ic J ( x - x Q ) 2 + (mh) 2 



where dP'/dx is the vorticity per unit length at the point : 
in the tunnel. 

.It will now" be assumed that the chord of the airfoil is 
reasonably small in comparison with the height of the tunnel. 
This being the case,. the approximation 



(x - x 0 ) 2 + (mh) 2 £T (mil)* 



is sufficiently precise for purposes of this analysis, and the 
term (x - x Q ) in the denominator of equation (38) may be 
neglected. 

The vertical velocity ( vt r )^ induced by all the images is 
then found by superposition as 

(y P )i = -L V Iliff / (x - x 0 ) ill dx 

jch 2 Z-x m 2 J . dx 

43=1 0 
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."ThP^'^f-^ (39) 

12h j dx 

o • 

This equation can be corrected for the effect of compress- 
ibility by means of expression (35). If, as was assumed, the 
chord of the airfoil is reasonably small as compared with the 
tunnel height, the polar angle cp of any point x 0 on the air- 
foil with respect to any point x on an image is nearly a right 
angle, so that in this case, the factor (35) is sensibly .equal to 



vfT- M 2 



The vertical velocity induced in a compressible stream by all 
the images is then 

c 

. v r « S — C (x - x 0 ) *£L ax (4-o) 

. 12h 3 ..\/ 1 - MV 
o 

t 

The vertical velocity v' b , induced at a point x p by the 
vortex sheet belonging to the airfoil itself, is given by equation 
(36) if T\ .and - v' b are substituted for P and v, respectively. 

The total vertical induced velocity v T at any point x Q on 
the airfoil in the tunnel is then the sum of v f r and v'^; that 
is, 



v = — * — / j — (x - x 0 ) 

2 * ~ Q [x - x Q 6h a (1 - u 2 y 



— dx. (41) 
dx 



The condition that the resultant of the induced velocity and 
the true axial velocity at the airfoil shall be tangential to the 
camber line requires that, at all points on ^the camber line, 

v j • dy c 

T « — - a' ■ _ _ (42) 
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where a r is the angle of attack of the airfoil in the tunnel; 
that is, the angle the chord line makes with the center line of 
the tunnel. The true velocity V rather than the apparent 
velocity V is used in equation (42), since the. vortex system 
used to represent the cambered airfoil in the tunnel is actually 
operating in a stream of velocity V when the airfoil thickness 
and wake are present, 

Relations "between Characteristics of Airfoil 

in Tunnel and in Free Air 

The preceding sections provide the basic information requir- 
ed for the development of relations between the characteristics - 
of the airfoil in the tunnel and in free air. The relations for 
the lift and moment coefficients and angle of attack are derived 
from the equations 1 of the preceding section by an extension of 
the method of Fourier series employed in Glauert's theory of thin 
airfoils (reference 12, pp., 87-93). To this end, the vorticity 
distributions for the airfoil in- the tunnel and in free air are 
each represented by a trigonometric series, the two series being 
similar in form but having undetermined coefficients. . By means 
of the equations of the preceding section, general relations are 
found between the coefficients of the two series. These general 
relations are then specialized to meet the requirement that the 
airfoil shall have the same value of the cotangent term of the 
series in the tunnel and in free air, this requirement being 
shown to be necessary to assure thatrthe essential characteristics 
of the pressure distribution will be sensibly the same in both 
cases/ By means of the relations between the coefficients ; 
expressions are then derived for the lift and moment coefficients 
and angle of attack of the airfoil in. free air in terms of the 
characteristics measured in the tunnel. The corresponding drag 
coefficient in free air can be* found from the drag measured in 
the tunnel by subtracting the pressure drag caused by the 
interference between the walls and the wake and referring the 
remaining drag to the true instead of the apparent dynamic 
pressure. Finally, a method is presented for correcting airfoil 
pressure distributions for the effect of tunnel walls. 

To carry out the analysis, points on- the airfoil are defined 
by a new coordinate 6 such that 

x = 1 c (1 - cos 0) (43) 
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and 



cbc = i c sin 0 d0 
2 



(44) 



The distribution of vorticity along the chord of the air- 
foil in free air is represented, after Glauert, by the trigono- 
metric series 



dx ! 



> A 0 cot - 9 + \ A n sin n9 
° 2 n 1 

n =i J 



(45) 



Equation (3*3) then gives the induced velocity at any point G on 
the airfoil as 



i =a/7- m 2 (-a 0 + y A n 



n=i 



cos nf) 



and equation (37) for the slope of the mean-camber line becomes 



dx 



= rc -a7 1 - M 2 A Q +a/~T~- M 2 A n cos nd 



n=i 



The coefficients are then given by the relations 

TT 



a -/vTl - M~ A Q = / —id da 

dx 



,2 



Air, 



A n 



V 1 _ M- ^ 



(46) 



(47) 



For the airfoil in free air the coefficients A n for n > 1 are 

thus functions of the camber-line shape only and are independent 
of the angle of attack. The coefficient A 0 is a function of 

both the camber-line shape and the angle of attack. 
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The chordvise lift distribution in free air is given by 



ft . P V *T = q 2 <LT 



ax 



dx 



V dx 



which after substitution from equation (45) can he written in 
coefficient form a3 



P = 



1 dL 
q dx 



r 

i -» 

4 < A 0 cot ± 0 + 

L 2 



1 

sin n0 J 



(43) 



Equation (48) illustrates the well-known fact that in free air 
the chordwise lift distribution consists essentially of two 
distinct parts. The one part, contributed by the sine terms and 
generally referred to as the basic lift distribution (reference 
13 ), depends in magnitude and form only upon the shape of the 
mean- camber line* The other part, defined, by the cotangent term 
and referred to as the additional lift distribution, is fixed in 
form and depends in magnitude up'on the angle of attack as well 
as upon the camber- line shape. 



The distribution of "vorticity for the airfoil in the tunnel 
is represented by 



JpLl = 2V \a q < cot i 0 +\ S A n T sin nG > 



dx 



L 



2 ^L. 



(49) 



J 



Substitution of this expression, together with expressions (43) 
and (44), into equation (41) gives, after integration, 

_ r /• * 



Ai 



1 - M 



2 



(2A 0 « + Ai') 



cos 9 



+ 1> 



1 cos n0 > 

J 
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where a is as defined by equation (2). Equation (42) for -the 
alopo of the mean- camber line thus becomes 



dy c , — 

-r-S = a' -N 1 

dx 



M' 



A 0 ' - 



M' 



As 



M 



2 ! A • 1 

I Ax' - 

L 



(2A 0 - 4 Ai') 



cos 0 



+\j 1 - M 2 A n ' cos n0 
The coefficients in this case are given by the relations 



17 



a' - *Jl - M ; 



0 i-M a l 0 2 As ; 



dx 



d8 



Tr 



Ai« - 2 (2A 0 « +, Ai ' ) = 1 1 I ilS. cos 6 a6 



1 - M' 



a/ 1 - M 2 "J ^ 



A n « * 



a/1 



M" 



TT 

r dy, 



-Uicos- n0 d0, n £ 2 
dx 



(50) 



>(S1) 



Thus for the airfoil in the tunnel the coefficient A 0 ? is a 
function of "both the angle of attack and the shape of the camber 
line , but the functional relationship is altered from what it was 
in free air by the inclusion of terms proportional to a. Further 
more, because of the appearance of the' term involving A 0 f .in the 
second of equations (51), the coefficient Ai 1 is in this case 
also a function of the angle of attack, as well as of the camber- 
line shape. Since A 0 ! appears in this equation multiplied by 
the factor a, the dependence upon the angle of attack is, how- 
ever, secondary as compared with the dependence upon the shape of 
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the camber line. As in the case of the airfoil in free air, the 
coefficients A n ! for n £ 2 are functions of the camber-line 
shape only. 

The chordwise lift distribution in the tunnel is given by 

ax,* ar» 2 dr y 

dx dx V dx 

or in coefficient form, 




In writing this equation the streamwise velocity gradient which 
results from the wall-wake interference (equation (24)) is 
ignored. It can be shown that the inclusion of this variable 
would give rise to correction terms of the order 003 T , Terms 
of this order are usually small as compared with the terms of 
order a. and t 03 ' considered in the theory and may therefore 
be neglected. 

It is apparent from equation (52) that, as in the case of 
the airfoil in free air, the lift distribution in the tunnel may 
be divided into two components. Now, however, the component 
which depends upon the angle of attach includes both the 
cotangent term and the first sine term. The component which is a 
function of the camber-line shape alone, comprises the sine terms 
corresponding to n £ 2. Again, these two components could be 
denoted by the terns "additional" and "basic" in the sense- 
previously employed; however, since the phrase "additional lift" 
already is so firmly established with reference to the distinctive 
cotangent term alone, this usage does not appear advisable in -the 
present case. For this reason, the terms of the series will be 
referred to by reference to their form or their position in the 
series, 

• Since it is the same airfoil which is being considered in 
both cases, equations (47) and (51) lead to the following general 
relations between the coefficients in free air and in the tunnels 
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a -\fl - M 2 A 0 = a' -tjl - Vf A 0 V + (k^ + - A 3 ' 

-\/l-M 3 ^ ' 2 



• A = A , - (2A 0 J .+ A 1 <) 

1 - M 



A 3 = V 



A n~ A n T 

In order to use these expressions to relate the charact er- 
istics of the; airfoil in free air with those in the tunnel,, it is 
necessary to choose some quantity or condition which will be main- 
tained the same in both cases and relate the remaining quantities 
in accordance with this choice. If it were possible, the- ideal, 
procedure would be to keep all the aerodynamic coefficients 
unaltered and to determine' a corresponding relationship between* 
the angle of attack in the tunnel and in f:ree air. To do this it 
would be necessary to keep all pressure and fractional forces the 
same in both cases, which can be accomplished only if the pres- 
sure distributions are identical. This would require that each 
of the coefficients in equation (48) be equal to the cor- 

responding coefficient A n * in equation (52). It is apparent 
from the second of equations (53), however, that this require- 
ment cannot, in general, be satisfied. 

Although the pressure distribution cannot be maintained 
completely unaltered in the transfer from the tunnel to free 
air, the general relations (53) can be specialized in such a 
way that the essential character of the distribution is un- 
changed. It is apparent that the component of lift contrib- 
uted by the first, or cotangent, term in equations (48) and 
(52) is different in form from that contributed by the series 
of sine terms. The cotangent component has an infinite value 
at the leading edge (8=0) and a relatively large chordwise 
gradient of lift over most of the chord of the airfoil. The. ■ 
sine-series component is finite at all points and, for air- 
foils ordinarily encountered in practice, has a relatively 
small chordwise gradient, except possibly in -the immediate 



r (53) 
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vicinity of the leading or trailing edges. The cotangent com- 
ponent with ita infinite pee.k pertains, of course, only to the 
hypothetical airfoil of infinitesimal thickness and zero leading 
edge- radius. For all real airfoils, the lift at the leading 
edge can never be infinite; however, even in this instance the 
lift distribution is characterized by a component the form of 
which is peaked near the leading edge and the magnitude of which 
varies markedly with the angle of attack. The magnitude of 
this component is a primary factor in determining the character 
of the pressure distribution, and even a relatively small change 
in magnitude may cause considerable change in the minimum pres- 
sure and in the chordwise pressure gradients attained on the 
surface of the airfoil. Further, the aerodynamic characteris- 
tics which depend upon these quantities, particularly the pro- 
file drag, maximum lift ; and critical compressibility speed, 
will be correspondingly altered. It follows that properly to 
correct airfoil data obtained in a wind tunnel to conditions 
in free air, the corrected quantities should correspond to the 
same magnitude of the, peaked lift component as exists on the 
airfoil in the tunnel. 

The requirement that the peaked component of lift on the 
real airfoil shall, be the- same in the tunnel , and in free air 
can be expressed with reference to the assumed airfoil of in- 
finitesimal thickness and camber by setting A Q equal to A Q ! 

in equations (53). The first of these equations, which re- 
lates the angle of attack in the two conditions, then -becomes 

a * a< + 2— (a 0 < + 1 A 2 ^ (54) 

VI - M V 2 I 

and the relations between the coefficients are 

A 0 « V 1. 

Al « Ai< (2A 0 ' +. A x ') 

1 - 

A 2 = A 3 < y (55) 



A - A 1 



t 
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Substitution from equations (55) into equation (48) gives 

P = 4 <A Q ' cot I 9 - -J^ (2A 0 » + AiO sinO + ^> A n * sin n9 I 
or L 



P = P» - (2A Q * + A^).sinS (56) 

Thus, if the angle of attack in the tunnel and the angle of 
attack in free air are 3uch as' to satisfy equation (54), the chord- 
wise lift distributions will differ by an amount defined by the 
second term on the right-hand side of equation (56). 

The lift, coefficient for the airfoil in free air is 



o 



which, after substitution from equations (44) and (48), can be 
integrated to give 

°l = *(2A 0 + A x ) (57) 

The quarter- chord-moment coefficient is 

l 

4 v ' 

o 

which becomes after integration 

J (Ai - A a ) (58) 



In usual wind-tunnel practice, the measured coefficients 
are referred to the apparent dynamic pressure q T . The lift 
distribution over the airfoil in the tunnel in terms of q ! is 

P i = JL _di^ _ 2 q. a P 

q' dx V q* dx 
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Substitution from equation (49) gives 



f ^ 
p» = 4 JL Ja 0 « cot i 9 + ^> A n " sin n9 ^ (59) 

q ^ n=i 



The lift and moment coefficients of the airfoil in the tunnel as 
referred to the apparent dynamic pressure are then, respectively, 



'I 

o 



= l P'd.(|) -« f," (2A 0 « +A x t) 



• <V = / P' d (^) - « (2Ao' + A x ') (60) 

and 



l 

vf * , (i-f) 4 (f)--^^ , -' A ' ,) (61) 

4 



4 

0 



Relations "between the coefficients in free air and in the 
tunnel can now be found -with the aid of equations (55). Sub- 
stitution of values from these equations into equation (57) 
gives 



Cj =„(2A 0 ' + A l .)(l-_H- pr ) 

4 q V 1 - M ' 

Substitution from equation (29) and neglect of correction terms 
of higher than the first order then give 

tfl * 2 I ( Mt ) 5 , A L2-(MQ 8 ]Cl+0>4(M<) a ] Tr , 

1 [ 1 " " [1-(M')^] 3/S " 1 - (MO 2 d 



c 2 = c z t, 1 



From equation (33) it can be seen that, to the first order 9 M 
may be replaced by M ! in this equation. The final equation for 
the correction of the measured lift coefficient is therefore 

o. = o, • (l - °— - f . > - ( M ') a . Ac 

1 1 [ 1 - (M' ) [T- (M') 2 J 3/2 

. r 2 . ( M»fj[i + o.4(MO^ ] (62) 

1 - (M') S 
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Similarly, substitution -of values from equations (55) into equa- 
tion (58) gives 



Ob = - i (Ax' - A 2 <) + * (2A Q - + A X ') c 



c 

~4 



1 - M" 



= 11 v t + 1 0 ■ _J2L_ 
a £ 4 1 - M* 



To the order of approximation previously employed, the final 
equation for the correction of the measured moment coefficient 
can be written 



r 

% = c m/ < X 



u c 

4 



£ 1 



Jli M j!_ A [2-(M') a ][l+0.4(M') 3 ] " 
[_l-(M')a]&/a iV0 " l-(M') 3 d 



+ V 



4[l-(M') 2 ] 



(63) 



The corresponding angle of attack in free air can be found from 
equation (54). Combination of equations (60) and (61) gives 



0 2 q I 2rt 



pZ* + *V""I 



L 



To the first order, equation (54) then gives for the corrected 
angle of attack in radian measure 



a = a' + 



cz' + ^mc'. 



(64a) 



or in degrees 



57.3 j 



r 



1 



a = a' + — ~~==. «; o,« + 4c * V 

at vi - (M>r L i J 



(64b) 



Numerical values of the compressibility factors appearing in 
equations (62) , (33), and (64) are given in table II. 
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It should not be implied from equations (64) that the general 
inclination of the stream at the position of the airfoil in the 
tunnel is actually different from what it would be if the walls 
were not present. The equations indicate rather that, with regard 
to the magnitude of the cotangent component of lift distribution, 
an airfoil at a given angle of attack in the tunnel behaves as 
though it were at a different angle in free. air. This difference 
occurs because the tunnel wails give rise effectively to a change 
in the curvature of the stream -at the position of the airfoil. 

As wa3 indicated previously, the essential character of the 
pressure distribution over a given airfoil will be the same in 
the tunnel and in free air, provided the magnitude" of the cotan- 
gent lift component is the 3ame in both cases; that is, provided 
the angles of attack 'are such as to satisfy equations (64), The 
exact shape of the pressure distributions, however, will still 
differ slightly for two reasons: (a) The interference between the 
lift and the tunnel walls causes a difference in chordwi3e lift 
distribution as required by equation (56), and (b) the. interference 
.between the wake and the walls gives rise to a longitudinal pres- . 
sure gradient defined by equation (25). The effect of these two 
influences upon the remaining airfoil characteristics } the profile- 
drag coefficient, must- be considered. 

As given by equation (56), the chord vise lift distributions 
in the tunnel and in free air differ by an amount 



AP = F* ~ P ■= ~~ ? ( 2Aq ' + Ai') sine 
1 - M 

which, by virtue of equation (60), may be written to the first 
order as 

AP = i _£ c « sin0 . (65) 

it. 1 - M 3 

The changes in peak pressure and pressure gradient brought about 
by this increment of lift distribution, unlike the changes which 
would accompany even a minor alteration of the . cotangent lift 
component, are ordinarily small. At low Mach numbers the drag 
depends primarily upon the character of the flow in the boundary 
layer, and, since this flow will not ordinarily be altered greatly 
by these small changes in the pressure distribution, the incre- 
ment of lift distribution should have only a small effect upon 'the 
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profile drag. At high Mach numbers the drag is determined prima- 
rily by the total-head losses in the shock waves which appear 
after the critical Mach number is reached; that is, after tne 
local speed of sound is obtained at the minimum pressure point on 
the airfoil. The critical Mach number is usually reduced by the 
change in peak pressure accompanying the change AP in lift 
distribution, but it can be shown that this reduction is ordi- 
narily very small. It is reasonable to expect that the change in 
profile drag at a given supercritical Mach number is correspond- 
ingly small. These changes are discussed in further detail later 
in the report, but for the present it may "be assumed that the. . 
difference in chordwise lift distribution between the tunnel and 
free air has only a negligible effect upon the profile drag. 

For usual airfoils and drag coefficients, the longitudinal 
pressure gradient defined by equation (25) is also small, and its 
effect upon the boundary- layer flow and hence upon the friction- 
al drag of the airfoil may be neglected. It will, however, in- 
crease the pressure drag by an amount which is comparable to 
differences already retained in the corrections to the lift and 
moment. Thi3 increase in pressure drag must be subtracted from 
the drag measured in the tunnel to obtain the true profile drag ' 
of the airfoil in free air. 

Glauert has shown (reference 2 , pp. 62-83) that in an 
iiicompressible fluid the drag experienced by an airfoil a3 the 
result of a streamwise pressure gradient is, in the notation of 
this paper, 

jt A ? dp 6h 2 .dp 

AD = Ac — = Ao — (66) 

3 dx jc dx 

In appendix A, it is shown that this relation is unchanged by the 
effect of fluid compressibility. Substitution of dp/dx from . 
equation (25) then gives for the drag due to the interference be- 
tween the wake and the walls 

AD = 04. ft« of 1 + OMWllj AO 

LQl - (M.) 2 J 3/ !f 

The time profile drag of the airfoil in free air is then 
D = D' - A.D 

-o a . ft' c [l - l + °- 4(M l )3 / Aaj 
d ^ 1 [l - (M') 2 ] 3/C 
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and the corresponding drag coefficient referred to the true 
dynamic pressure is 

d . *c d WV [1 - (M') 2 ] 3/H J 

Substitution from equation (29) give3 for the final correction to 
the measured drag coefficient 



Cj = C j ' 



1 . 2 - (M 'T Aa- 1 + °- 4 ^'r Aa 
[l - (M') s 1 3/a [l - .(M-)3p/2 

_ & - (M') 2 ][l + 0.4(MQ»] t ^ 1 (6?) 



1 - (M«) 



'J 



It i/iil be noted that, of the two correction terms involving A a 
in this equation, the first appears as a result of the change in 
dynamic pressure occasioned by the interference between the walls 
and the airfoil thickness; the second represents the effect of 
the pressure gradient induced by the interference between the 
walls and the wake. The correction term containing .t c^ 1 appears 
as a result of the change in dynamic pressure caused by the wall- 
wake interference. Numerical values of the* functions of M* 
which appear in equation (67) are given in table II. The cor- 
rected drag coefficient corresponds, of course, to the corrected 
lift and moment coefficients as given by equations (62) and (63) 
and to the corrected angle of attack as given by equation (64a) 
or (64b). 

The drag correction of equation (67) was determined partic- 
ularly for drags measured with a balance and, as derived, is not 
necessarily correct for drags measured by the wake-survey method. 
It can be shown, however, from theoretical considerations of 
momentum and continuity in a two-dimensional-flow tunnel that for 
normal ratios of airfoil chord, to tunnel height , the ordinary 
type of wake survey derived for free-air conditions gives, when 
applied in the tunnel, a value of the drag equal to that measured 
by the balance except for a negligible difference of less than 
one-half of one percent. Equation (67) may thus also be used to 
correct drag coefficients determined. by the wake- survey method. 

It should be noted that no correction to the drag has been 
made for any pressure gradient which may exist inherently in the 
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tunnel as- a result of the streamwise growth of the boundary layer 
on the tunnel walla. Most modern tunnels are constructed so that 
this pressure gradient is sensibly zero; however, if such a gra- * 
dient does exist and its magnitude is known, an approximate cor- 
rection to the airfoil drag can be made by means of equation (66). 

There remains the necessity for correcting the measured' pres- 
sure distribution over the surface of the airfoil. The' pressure 
at any point on the airfoil is conveniently expressed by the pres- 
sure coefficient Si ' defined by 

H - p 7 
Si -± 

or by the pressure coefficient P z> defined by 

p, - Vl : p (69) 

1 i ■ 

where p^ is the local static pressure *on the surface of the air- 
foil and H, p, and q are, respectively, the total head, static 
pressure, and dynamic pressure of the undisturbed stream. As 
indicated in reference 14, in a compressible stream, 

H = p + q (1 + tl) (70) 

where (1 + r|). for air (7 = 1.4) is defined by the series 

1 + Tl ss 1 + + + + . . . (71) 

4 40 1600 

M being the Mach number of the stream. From these relations it 
is readily shown that 

S z = (1 + Tj ) - P z (72) 

A curve of (1 -f i}) versus M, as calculated from equation (71), 
is given in figure 4. 

In reference 6 a method is presented for the determination 
of the pressure distribution around an airfoil in an incompress- 
ible stream when the lift distribution along the chord and the 



(68) 
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pressure distribution over the base profile are known. The upper- 
and lower- surface pressures at any chordwise station x are given 
in coefficient form by 



= 1 



1 

Pi, - P 



(i - P f ) 

[d-P f )-^p] 2 
(i - Pf ) 



2 



(73) 



where P^ is the pressure coefficient on the base profile at x, 

and P is the coefficient of lift per unit of chord at . x. By 
following the basic reasoning of reference 6 and assuming that 
the induced velocities at the surface of the airfoil are small as 
compared with the velocity of the undisturbed stream, it is readily 
shown that equations (73) may also he applied to the pressure 
distribution in a compressible stream. In such application, the 
values of Tjj, P*l, Pf, and P must all correspond, of course, to 
the same free- stream Mach number. 

The measured pressure distribution is now readily corrected 
for the effect of the tunnel walls. It i3 only necessary to refer 
the measured pressure coefficients to the true instead of the 
apparent dynamic pressure and remove the effect of the lift dis- 
tribution represented by equation (65). Strictly speaking, cor- 
rection should also be made for the pressure gradient due to the 
wall-wake interference; however, in practical. tests such cor- 
rection is small and may be neglected. The detailed procedure is 
then as follows: 

(l) The apparent upper- and lower-surface pressure coef- 
ficients 

H'- p L rV ' H - p L » 

= Si- and S T 1 = — 

u q' L q» 



are obtained from the experimental results for the various chord- 
wise stations. 
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(2) These pressure coefficients are referred to the- true 
dynamic pressure by means of the equations 



S L *=S T ■ fl - 2 -< M ' )2 ,Aq. [2-(M') 2 jLU0.4(MQl ^ 

L L vJ L 1 [i-(m')T /0 - . 1 i-^') 3 



> (74) 



(3) The quantities (1 - Pu*) and (1 -Pi*) are determined 
in accordance with equation (72) as 

1 - P n * = Sjf - r, . . , 

(75) 

1 . r L * = s L * - n 

where ti is determined by figure 4 for the true Mach number as 
given by equation (33). 

(4) The chordvise lift distribution in the tunnel i3 found 

from .' 

p* = SLlJSL = Su * - S L * (76) 
q 

(5) The chordwise lift distribution in free air is determin- 
ed from equation (65), which may be written 

P = P* - — CT , - a ?e «Z' (?7) 



1 - (M«) 



vhere P e is given by 

4 



P e = - sin 6 
v ' it 



This quantity 9 vhich is termed the. " interference lift distribu- 
tion/' is seen to be elliptic in form. Values of ' P.^ at stand- 
ard chordvise stations are given in table III . 
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(6) The quantity (1 - Pf), 'where Pf is the base-profile 
pressure coefficient in free air, is given "by the equation 

' :' (i . Pf) •. a . Vr./jgZg *_ sgZg 

which is obtained by combing equations (73). . 

(7) The values of P and (1 - Pf) being known, the upper- 
and lower- surface pressure coefficients Py -and P^ are deter- 
mined from equations (73)/ If. desired, the corresponding coef- 
ficients Su and. Sl can be found from equation (72), 

The corrected pressure distribution obtained, by this method 
corresponds to the corrected angle of attack as given by equation 
(64a) or (64b) and to the corrected, lift and moment coefficients 
as given by equations (62) and (63), 

It has been mentioned previously that the correction to the 
angle of attack appearing in equations (64) does not represent an 
actual rotation of the stream direction. This fact is implicit 
in the derivation of the equations, but it can also be demon- 
strated 'by simple considerations of force and momentum. For this 
purpose it is sufficient to consider a simple incompressible 
potential flow in the. tunnel and ignore the effect of the profile 
drag. Assume for the time being that, because of the interference 
between the airfoil and the tunnel walls, the general direction 
.of the. stream at the airfoil- is inclined from its original 
direction parallel to the tunnel walls. For potential flow the 
resultant force acting on the airfoil must be at right angles to 
the local direction of the stream. The airfoil thus would be 
acted upon under the assumed conditions by a component of force 
parallel to the center line of the tunnel and would in reaction 
exert an equal and opposite force on the floy. Since the. tunnel 
walls cannot in a potential flow exert a force parallel to the 
center line, this - longitudinal force would have to be balanced 
by a difference of pressure or momentum between two stations in 
the tunnel, one upstream and one downstream from the airfoil. If 
the stations are taken far enough from the airfoil that its ' 
induced velocities are negligible, conditions across, the funnel " 
are uniform at each station. It then follows, from considerations 
of continuity of the incompressible flow. in the tunnel that the 
conditions at the two stations are identical, and no difference 
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of pressure or momentum is possible. Thus the original assump- 
tion of a general rotation of the stream direction at the position 
of the airfoil is untenable. This conclusion is not changed by 
the effects of fluid compressibility. Furthermore, the fact that 
the introduction of the profile drag and the accompanying wake 
causes a pressure difference between the two stations likewise 
does not alter the result, as the wake effects are considered in 
the theory to be superposed on the potential- flow field. Thus, 
the angle correction appearing in equations (64) must be due to 
some cause other than a general inclination of the stream. As 
previously pointed out, it is actually due to an effective change 
in the curvature of the stream at the position of the airfoil and 
is a direct consequence of the requirement that the airfoil in 
this stream shall have the same cotangent component of lift 
distribution as does the airfoil in free air. These considerations 
are important in the proper interpretation of drag measurements 
from a two~dimensional-f low" tunnel. 

In the development of the correction to the measured drag 
coefficient, it was assumed that the increment AP in chordwise 
lift distribution between the tunnel and- free air has only a 
negligible effect upon the profile drag. A better idea of the 
nature of the effect .can be had by further examination of the 
difference between the two cases. It follows from equations (55) 
that, if the angles of attack in the tunnel and in free air are 
related as required by equation (54) or (64), the transposition 
of a given airfoil from free air to the tunnel is equivalent to 
increasing the coefficient A x for the airfoil in free air by 
an amount 



which can be written to the first order as 



As can be seen from equation (46), this can be accomplished by 
maintaining the angle of attack unaltered in free air and chang- 
ing the ordinate of the mean-camber line at every point by an 
amount Ay c such that 



AA, = Ai' - Ax = ° - (2A Q < + Ax') 
1 - M 



d(Ay c ) 




1 



a 



(80) 
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The value of Ay as a fraction of the chord is : then 

d(Ay c ) 




dx 



00 



+ C 



which after substitution from equations (44) and (80) can be inte- 
grated to obtain 



Ay c _ 1 a 



ci cos 2d + C 



The constant of integration " C is determined by the condition 
that Ay c /c = 0 at 9 = 0 'and 9 = it.. The equation for the 
change In the camber line then "becomes finally 



- 1 * C7 



(f)-(iTj 



(31) 



.This is the equation of a parabola with vertex at* the midchord 
point and has the same form as the equation for the camber line 
of an MCA conventional airfoil with maximum camber at the mid- 
chord point (reference 15). The maximum change in camber is 



= J-_EL_c- 7 (82) 
V c 4ax 4 * \ITZ~p 

Thus, if the angles of attack of the airfoil in the tunnel and 
free air are' adjusted as required by equations (64), the wall 
interference in the tunnel has the same effect upon the chord- 
wise lift distribution as would an increase in camber in free air. 

As a possible instance of a test for the determination of 
the drag of an airfoil of large chord at a low Mach number and 
low lift coefficient, consider the case of an airfoil in a tunnel 
providing a chord-height ratio of 0.5. The value of a. is /then 
6.051. Assume that the angle of attack a 1 in the tunnel is 
adjusted as required by equations (64) to correspond to an angle 
a giving a lift coefficient c^ of 0.3 in free air. Assuming 

that the Mach number is sufficiently low that the effect of 
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compressibility may be neglected in computing the tunnel-wall 
corrections, the change of maximum camber required in free air to 
duplicate the effect of the tunnel walls is given by equation (82) 
as 



(^) ,= JL (0.051)(0.3) = 0.0012 
^ c /max' 4,c 

An estimate based upon experimental data has been made of the 
effect upon the profile drag of a change in camber of this magni- 
tude for" an NACA conventional airfoil of moderate camber and 15- 
percent thickness with maximum camber at the midchord point. The 
result indicates that neglecting the effect upon the profile drag 
of the change in lift distribution caused by the tunnel walls 
leads in this case to an error in the final corrected drag co- 
efficient of less than 0.0001. This is within the usual limits 
of experimental accuracy. The correction terms included in equa- 
tion (67) amount in this instance to approximately 0.0004. If 
the chord-height ratio were increased to 1.0, the error in the 
drag coefficient would be increased to 0.0004, which is well out- 
side the limits of 'experimental accuracy. This indicates the 
desirability of limiting the chord-height ratio if accurate 
measurements of the profile drag are desired, even at low values 
of the lift coefficient and Mach number/ At higher values^ of the 
lift coefficient or Mach number the permissible chord-height 
ratio must be reduced correspondingly. 

The foregoing comparison is based upon, the specific case of 
an airfoil with maximum camber originally at the midchord point 
and is not necessarily applicable to other typos of airfoils. 
For families of airfoils which have a smaller variation of drag 
with camber than do the MCA conventional sections, the error 
introduced by neglecting the effect of the change ,ih lift dis- 
tribution is correspondingly less. In any event, satisfactory 
accuracy can be obtained in the measurement of drag at low' lift 
coefficients and Mach numbers by keeping the chord-height ratio 
within a suitable limit - say 0.7. A possible exception is an 
airfoil having an essentially flat pressure distribution in the 
region of transition from laminar to turbulent flow in the bound- 
ary layer. In such a case the changes in pressure gradient may 
shift the point of transition and noticeably alter the profile 
drag; however, for any such sensitive. airfoil, alterations from 
this source are no more serious than similar changes which may 
accompany the small variations in pressure distribution caused in 
any practical application by irregularities in construction. 
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Some measure of the .effect of the increment • AP in chord- 
vise lift distribution upon measurements of airfoil characteristics 
at high Mach numbers can be obtained by calculating the .change in 
critical Mach number caused by this increment. Such a calcula- 
tion has been made for an airfoil with minimum pressure originally 
at the midchord point. Since the increment AP is a maximum at 
midchord, this represents the worst possible case, as regards the 
change in critical Mach number.' For a chord -height ratio. of 0.25 ; 
which is considerably larger than that ordinarily employed in 
tests at high Mach numbers , the critical 'Mach number was fotind to 
be reduced by approximately 0.00.1 at a lift coefficient of 0.3. 
A change of this magnitude is insignificant. It may be expected 
that the accompanying change in the aerodynamic coefficients in 
the vicinity of the critical Mach' number will be correspondingly 
small, 

/ 

THE PHENOMENON OF CHOKING 

Consider the compressible adiabatic flow of a fluid in an 
elementary stream, tube of varying area Aj, as shown in figure 

5(a). Continuity of flow requires that the product PjVjAj be 
.constant, where p^, Vj, and Aj are the local values of- 
density, velocity, and area, respectively > at any station along 
the tube. In consequence, the logarithmic derivative must vanish; 
that is, 



dp 7 dV 7 dA 7 

JJL + —I + _J = 0 (83) 

p z v z h 

Bernoulli's equation for compressible flow requires that 

• 5l» ~V z dV 2 (84) 

where p^ , is the local pressure. Defining V c as the local 
velocity of sound, then 
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30 equation (84) becomes , after substituting the value of dpj 
in that equation, 

vhere Mj is the local Mach number. ' 

Substituting this relation into equation (83) gives 

• ' dV, dA 7 

(1 , M a ) __±. = 1 (85) 

1 V 2 . A l 

From this well- known relation it is seen that' at subsonic 
speeds the usual "behavior associated with incompressible flow is 
obtained; namely, that as the area increases the velocity de- 
creases. At supersonic speeds, however, the behavior is reversed 
in that as the area increases the velocity increases. When the 
local Mach number is unity it is seen that dA = 0; that is, if 
the velocity of sound is attained in the tube it can only be 
attained where the area has its minimum value. 

When the local velocity of sound is attained at the minimum 
area section, the local Mach number at any other section, deter- 
mined by the ratio of the area at that section to the minimum area 
may be less or, in some cases, greater than unity depending upon 
the conditions promoting the. flow in the tube. The. nature of such 
flows can be studied by considering the change in the character 
of flow iri the stream tube of figure 5(a) as the downstream pres- 
sure p s is decreased with respect to' -the ' upstream pressure p 1# 
If pi - P 2 is small so that completely subsonic flow is main- 
tained' in the tube, the nature of the velocity variation along 
the tube is that usually associated with incompressible flow as 
seen in curve I of figure 5(b)-. When pi -. pa is' increased so 
that sonic speed is just reached in the. minimum area section, the 
variation of velocity along the tube -becomes that shown in curve 
II. Any further decrease of the pressure p 2 cannot alter the 
flow upstream of ' the minimum area, since the velocity at the 
minimum section cannot exceed the velocity. of sound. The only 
effect of decreasing the downstream pressure is to promote a • 
supersonic flow region downstream of the minimum area, as shown 
by curve III. of figure 5(b). This region is terminated by an 
abrupt return, through a compression shock; wave, to subsonic flow. 
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The position of this terminal shock wave must, be such as to bring 
about the necessary conversion of kinetic to thermal . energy that 
is required to promote the downstream pressure ps. For present 
purposes, the most important point concerning the flow as describ- 
ed is that when the velocity of sound is attained' at the minimum 
area section, no further increase in the flow rate can be, made 
regardless" of the extent of the supersonic flow region downstream 
of ' this section. When this maximum flow rate has been reached, 
the stream tube is said to be "choked." 

What has been said concerning the choking of a single stream 
tube applies to the complete system of stream .tubes past an air- 
foil mounted in a two-dimensional-flow tunnel, a3 shown in figure 
6. That is to say, when the* velocity of the undisturbed flow far v 
upstream in the tunnel reaches a certain value, sonic velocity is 
attained at the point of minimum area of. each elementary stream 
tube between the airfoil and the upper wall of the tunnel. It is 
important to note that the locus of the points of. minimum area of 
the separate stream tubes does not necessarily coincide with the 
shortest line between the airfoil and the upper wall. This is 
illustrated in figure 6, where the line A represents , the short- 
est distance between the airfoil and the wall. If the conditions 
of flow were uniform across the stream at each chordwise station, 
the flow between the airfoil and the wall would be the seine as in 
a .single elementary stream tube, and sonic velocity would neces- 
sarily, be attained along line A, In the actual case, however, 
the flow is two-dimensional, and sonic velocity is attained along 
some line, such as line B, not coincident with A, A similar 
situation exists in the space between the airfoil and the lower 
wall of the tunnel, where the sonic velocity. is attained along 
some line D. . As before, this line does not necessarily coincide 
with line C, the shortest line which can be drawn from the 
lower surface of the airfoil to the lower wall, (In order to 
avoid an apparent contradiction with the requirements of continu- 
ity, it must be kept in mind that the velocity vector, is not, in 
general, perpendicular to either lines A- and C* or B and C ) Sonic 
speed is generally not attained coincident ally along lines B and 
D, Once it is attained along both these lines, however, the rate 
of flow past the airfoil in the tunnel can undergo no further 
increase. The Mach number of the flow ahead of the. airfoil then . 
has its maximum attainable value. This value is described as the 
"apparent choking Mach number. ! ! 

In practice, the lines of sonic speed lie very close to the 
lines defining the shortest distance between the airfoil and the 
tunnel walls. For purposes of analysis, it will be assumed that 
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they are coincident, that is, that lines' E and D coincide, respec- 
tively, with lines A and C, Under these conditions, the calculated 
rate of flow in the tunnel (which mu3t in any event be equal to 
, the rate of flow across lines A and C) will be somewhat greater 
' than that which actually exists when the lines B and D have their 
true positions. The assumption of unidimensional flow will thus 
lead to a computed choking Mach number, which is slightly greater 
than the theoretically correct value. 

On the basis of the foregoing assumption a relationship be- 
tween the model size and the choking Mach 'number can be obtained 
from elementary considerations. The velocity V 1 and density 
p T of the flow far forward of the model, where the cross- 
sectional area is A f , are constant across the stream. The 
velocity and density p m across the sonic- speed lines B and 

D of figure S, where the area has the minimum value A m , are again 

constant across the stream. The velocity V m is the : local sonic 

speed V c ^ so that the equation of continuity becomes 

p'VA' - PmV^Aa 

Assuming adiabatic relations, the density and velocity terms 
can "be related to the Mach number far upstream, which is now the 
apparent choking Mach number. The end result is that the ratio 
of the area of the undisturbed stream to the minimum flow area 
can he expressed in terms of the apparent choking Mach number 



M 'ch as 



7+1 



Am M 'ch 



+ m j_ (M ' ch) " 1 ^ 



(86) 

J J 



The area ratio is clearly 



h 



A m . h - t. 



'P 



where h is the tunnel height and t p the projected thickness 
of the airfoil normal to the flow direction. For reasons which 
will he evident later, the projected thickness in this relation 
will he replaced by an "effective" thickness t e . 
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Taking the value for y for air as 1.4, equation (86) be- 
comes 



In figure 7, the ratio t e /h is plotted as a function of 

the apparent choking Mach number. The region above the curve 
represents an impossible state of flow, As a matter of interest 
the results are shown for the supersonic- as well as the subsonic- 
flow regime, although for the purpose of this report only the 
subsonic choking Mach numbers will be considered. 

In writing equation (87), the projected thickness was re- 
placed by an effective thickness. If choking occurred as was 
assumed in the preceding analysis, then the effective thickness 
determining choking would be, 'of course, the projected thickness. 
In any real case, although the. effective thickness may never be 
less than the projected thickness, it may be greater for two 
reasons. First, if the angle of attack is sufficiently large in 
absolute value, one of the lines B or D.may move downstream of 
the trailing edge because of the continued contraction aft of the 
airfoil of the portion of the stream passing that line. Second, 
since on any aerodynamic body there exists, because of the action 
of viscosity, a boundary layer wherein the velocity must be re- 
duced below the velocity in the otherwise unaffected flow field, 
it follows that the velocity of sound cannot be attained at those 
points close to the airfoil surface on the lines B and D of 
figure 6. 

To estimate the choking Mach number in any practical case, 
it is necessary to assume that the effective thickness is equal 
to the projected thickness of the airfoil. Because of the 
possible contraction of part of the stream aft of the airfoil, as 
well as of the assumption. that unidimensional flow exists as 
previously described, this procedure will lead to a computed 
choking Mach number which is greater than the theoretically 
correct value for an ideal, incompressible fluid. Further, the 
influence of the boundary - layer will cause the actual choking 
Mach number to be even less than this theoretically correct value. 
Thus the use of the projected thickness in the computation may be 
expected to lead to an overestimation of the choking Mach number. 




(M' ch ) a - 11 3 



(87) 
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' T'ho effect of the boundary layer in this regard may beet be 
illustrated by the .case of a flat plate set at zero angle of 
attack in a two-d3men3ional-f low- wind tunnel. Since the project- 
ed thickness is zero,, the previously developed theory would indi- 
cate that no choking would occur in thia case. Actually, because 
of the fact that the plate has a boundary layer and an accompany- 
ing wake, choking does occur, as is shown in the following dis- 
cussion, . . 

It yc.s seen in the section on wake effect, wherein the effect 
of confining the wake of 1 a body experiencing drag was considered, 
that when the influence of. the wake spreads to the walls so that 
a vavlform velocity field again exists, the temperature at this 
downstream position is related to the temperature upstream of the 
model by equation (13), using equation (9), the ratio of the 
corresponding velocities may be seen to be 



V". 



l + 7(K') a fl -^fcY| 
L 2 \h/J 



V" (7 + 1){M ? ) 



\2 



. " " (7 + l)(M') a ' ' 

The velocity ratio is* imaginary when the sign of the group 
of terms under the radical is negative. The functional relation- 
ship between the choked Mach number r.nd the drag- density factor 
Tc d 5 , found by equating the terms under the radical to zero and 
solving the resulting equation, is thus determined as 

- ■ 1 + 7&'<af i /r'r 1 "^'? ! (69) 

T cd T = * r— • < 1 - a /I - i , r~< f l yy >> 

.27(M. ch ) 2 ^ |_1 + 7(&' ch )*j J" 

where, as before. T = ~ / 2 Y, Setting 7 = 1.4 for air gives 

. 4 \h J 



d 2.8(M' oh )« | N 



(M»ch) 



l+1.4(M' cll ) 2 



(90) 



a graph of ' this function is shown on f igure 8, The effect of drag 
on choking for supersonic as well as subsonic wind tunnels is 
shown as a matter of interest. 
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The manner in which drag promotes choking may be comprehended, 
by examining the variations of the ratios V"/V r ' and T'y'T* in 
equations (88) and (13) as the value of Tc d ' is increased. In 
the case of the subsonic wind tunnel, the effect of increasing 
Tea 1 is to increase V'/V, On the other hand, T ,f /T f and hence. 
V c "/ V c ? are reduced. Consequently, M"/M* is increased. In 
the case of the supersonic wind tunnel the .effect of increasing 
7cd ? is to decrease V'/V 1 and to increase T"/T x and hence 
V 0 n / v c I j Consequently, M"/M* is reduced in this case. In both 
cases choking occurs when the value of T c^ 1 is such as to make 
the downstream Mac'h maaber M" unity. 

There is one definite limitation of the -previous analysis in 
that it was assumed that the effective tunnel area remained con- 
stant at least until the wake had spread to the walls so that uni- 
form flow conditions were obtained across the stream. Such a 
condition does not prevail in any conventional wind tunnel, never- 
theless the results are useful in providing approximate values for 
the effect of drag as it determines choking. For example, a flat 
plate having an apparent drag coefficient of 0.007, if the chord- 
height ratio were 0.5, would choke a subsonic wind tunnel at a 
Mach number of 0.95 df choking occurred as assumed* in the analysis. 
The serious influence of drag on choking for airfoils for which 
the drag coefficient may be many times this value is evident. 

To summarize, it has ( been shown that choking will occur in a 
wind tunnel' as a result of the confinement of the flow caused by 
the presence of the model and its wake. In the case of airfoils 
: of normal thickness, choking will usually be determined by. the 
effective dimensions of the body - that is, by the actual dimen- 
sions modified for the effects of boundary layer and stream con- 
traction aft of the airfoil as previously described. Properly, 
the boundary- layer effect is a. drag influence, but since its 
contribution i3 usually small it is most convenient to classify 
such confinement effects along with those due to the physical air- 
foil dimensions. In the case of very thin airfoils, at small • 
angles of attack, choking will usually result from the confining 
effect of the wake rather than the effect. of the airfoil thickness. 

Once the choking Mach number is reached, no further increase 
in tunnel power can affect the apparent Mach number. Such an in- 
crease will only serve to extend the supersonic flow region down- 
stream of the lines of sonic speed. The forces experienced by 
the airfoil at choking thus vary depending on the power input to 
the wind tunnel ^ 
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As a final consideration it should be noted that tho flow in 
the tunnel at choking does not correspond to any real flow over 
an airfoil in free air, " Since the choking Mach number approaches 
unity as the tunnel height h "becomes infinite, flow in the 
tunnel at choking, if it is to correspond to any flow in free air, 
must correspond to the flow that would occur , around an airfoil in 
a free stream, moving at the velocity of sound. It can 'be demon- 
strated, however, % that such a correfjpondeii.ee is impossible. 
Experimental evidence indicates that the flow conditions existing 
in tho tunnel at choking are essentially steady state. That the 
flow about an airfoil in a free stream having the velocity of 
sound cannot be a steady- state flow can be readily shown, For 
instance, it was demonstrated previously that in any stream tube 
the velocity of sound, if it is attained at all, must be attained 
at the minimum area section. That is to say, the rate of flow 
per unit area is a maximum, whore the velocity is the velocity of 
sound. Now, presuppose a steady-state flow in the stream tubes 
in the vicinity of an airfoil in free air when the stream 
velocity is sonic speed. If tho velocity either increases or 
decreases as the flow passes the airfoil, the stream tubes mu$t , 
expand. This is clearly impossible, since" the disturbance to the 
flow would then increase continuously as the distance from the 
airfoil increases. On the other hand, if the velocity remains 
the velocity of sound in each stream, tube, the streamlines will 
then have the same shape at all distances from the airfoil. Also, 
the pressure will remain constant throughout the entire flow 
field. This is, of course., .impossible, since pressure differences 
are necessary to promote the required changes in the direction of 
flow past the airfoil. A steady-state flow similar to that observ- 
ed in the tunnel at choking therefore cannot exist in free air at 
a free- stream Mach number. of unity, Thus at the choking Mach 
number, the flow at the airfoil in the tunnel cannot correspond 
to any' flow in free air. It follows that, at choking, tho 
' influence of the tunnel walls cannot, be corrected for. Further, 
in the range of Mach numbers close to choking, where the flow is 
influenced to any "extent by the incipient choking restriction, 
any correction for wail interference must be of doubtful validity. 

That the flow at or close to choking cannot be corrected for 
the interference effects of the tunnel walls may be reasoned from 
another, point of view. The assumption that it is permissible to 
correct wind-tunnel test data for the influence of the walls is 
justified only when the influence on the flow near the model is 
of such a uniform nature' as not to alter the general character of 
the flow materially from some corresponding flow, in free air. 
For instance, a velocity correction for wall interference may be 
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applied with confidence only if the velocity increment resulting 
from such interference is constant or nearly constant over that 
portion of the flow field wherein the influence of the model on 
the flow is important. Viewed in this light, it is clear that 
at or clone to choking no correction can properly "be applied, 
since an Important influence of the model on the flow is felt 
over a range extending close or up to the walls, within which 
range the influence of the walls on the flow is not at all uniform. 

It is thus clear that the equations which have been derived 
for correcting the. test data obtained in a subsonic two- 
dimensio?:±al-fIov wind tunnel for 'the effects of wall interference, 
cannot apply at the choking Ma oh number nor for a range of Mach 
numbers below the choking value. Moreover, when the model is not 
symmetrically disposed, the flow will, in general, attain sonic 
velocity across the stream on one side of the airfoil. before it 
does on the other. In such cases, it is' to be expected that the 
range of Mach numbers below choking for which the corrections are 
invalid is extended over that which would occur with a more 
nearly symmetrical flow pattern. 



DISCUSSION 

There is, at present, only a very limited amount of ex- 
perimental data available which can be considered satisfactory 
for determining the accuracy of the theoretical interference 
corrections derived in this report. Moreover, none of the 
available data were obtained at sufficiently high Mach numbers 
to permit an evaluation of the accuracy of the theory with regard 
to the effect of compressibility. 

- In figure 9 are shown the experimentally determined varia- 
tions of lift coefficient with angle of attack for several NACA 
0012 airfoils, having different chord -Iteight ratios. The data 
foi" those models for which the chord-height ratios are 0.25, 
0.5, and 0.8 were obtained from tests in the 7- by 10-foot wind 
tunnel at the Ames Aeronautical- Laboratory. These models were 
of 6-foot span mounted across the 7-foot dimension of the test 
section; 6- inch- span dummy ends were used in an attempt to obtain 
two-dimensional flew. A gap of about 1/32 inch occurred between 
the test panel , which was connected to the balance frame, and the 
dummy ends, which were fastened to the tunnel walls. The lift was 
determined both from force tests and by integration of chordwise 
pressure distributions at a section close to midspan. The data 
presented here are those obtained from the pressure distributions. 
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The data for. the model for which the chord -height ratio is 1-0' 
were obtained from teste in the low- turbulence wind tunnel of the 
Langley Memorial Aeronautical Laboratory. This 3-foot- span model 
was fastened directly to the side walls of the tunnel such that 
no air gap existed , and the lift was determined from measurements 
of the reaction on the roof and floor of the test section. The 
test results for zhe various models are shown, uncorrected for 
tunnel-wall Interference, in figure 5(a), In figure 5(b), are 
shown the same data corrected for wail interference by means of 
equations (62) and (34b). For all the .models, the correction 
term depending upon Tc^' is negligibly small. The test 
Reynolds numbers range from 2,000,000 to 6,000,000. It is seen 
that the corrected data obtained with the models for which 
(c/h) equals 0.25, 0.5, and 1.0 agree well with one another and 
with the 'section lift characteristics as obtained from tests in 
the MCA variable-density wind, tunnel (reference- 15). The data 
obtained with the model for which (c/h) equals 0.8, when 
corrected 5 indicate a lower lift- curve slope than do the other 
data. This is thought to be clue to the effect of air leakage 
through the gaps at the ends' of the test panel, the influence of 
which may be expected to become more pronounced as the chord of 
the airfoil is increased relative to the span. 

In this regard, unreported tests in the Langley low- 
turbulence wind tunnel have shown that the presence of any gap 
through which leakage can occur will influence .the aerodynamic 
characteristics to a surprisingly marked extent. This fact was 
also demonstrated by the Ames Laboratory tests on the NACA 0012 
airfoils. A comparison of the lift characteristics obtained from 
balance measurements with those derived by integration -of the 
pressure distributions, which are those given in figure 5, showed 
the lift- curve slopes for the former to be definitely lower than 
those for the latter, This indicates that 'the lift near the 
center of the test panel exceeded that at the sections near the 
gaps; that is, that the flow was definitely not two-dimensional. 

In figure 6(a) is shown the experimental variation of lift 
coefficient with angle of attack for an MCA 23012 for which 
(c/h) equals 1.0. These data were obtained in the Langley low- 
turbulence wind tunnel at test Reynolds numbers of 4,560,000 and. 
6,450,000. The same data corrected for tunnel-wall interference 
by means of equations (62) and. (64b) are shown in figure 6(b), 
together with section lift characteristics as obtained in the 
variable -density wind tunnel at an effective Reynolds number of, 
5,000,000 (reference 16). The corrected data are seen to be in 
excellent agreement with the results from the variables-density* 
tunnel . 
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In figure 7(a) is shown the variation of quarter- chord - 
moment coefficient with lift coefficient for the NACA 0012 air- 
foils as obtained from the 7- by 10-foot vind- tunnel tests 
previously described. In figure 7(b) are shown the same data as 
corrected for the interference of the tunnel walls "by means of 
equation (62) and (63). The section moment characteristics for 
this airfoil as obtained from tests in the variable-density wind 
tunnel (reference 15) are also shown for comparison. It is seen 
that the corrected data are in fair agreement with the data from 
the variable-density wind tunnel, except for the model for which 
(c/h) equals 0.8, It is believed that this disagreement is 
again due oo the effects of air leakage through the gaps between 
the test panel and the dummy ends, and not to any shortcoming in 
the theory. 

In figures 13 and 13, the uncorrected and corrected profile- 
drag coefficients for six symmetrical bodies at zero angle of 
attack are plotted as a function of the experimental chord-height 
ratio. The uncorrected experimental values 03 s , shown by the 
crosses , are taken from results reported by Fage in reference 17. 
The theoretically corrected values ca, indicated by the circled 
points, were computed from equation (67) for M 1 = 0. The extra- 
polated free-air value given in reference 17 for each of the 
bodies is indicated by a horizontal dashed line. It is seen that 
the corrected points are in good agreement with the extrapolated 
free-air values. In view of : the assumptions made in the theoreti- 
cal development/ the relative accuracy of the corrections at 
large chord-height ratios and large drag coefficients is remark- 
able, particularly in the case of the circular cylinder. 

Glauert (reference 2, pp. 56-57) suggests for the drag cor- 
rection in an incompressible fluid a formula which may be written 
in the -notation of this paper as .< 



where (t/'c) is the thickness ratio of the airfoil. In this 
equation, as in equation (67), the first correction term appears 
as a result of the interference between the airfoil thickness 
and the tunnel walls, and' is identical with the corresponding 
term in equation (67) for M 1 = 0. The remaining term is an 
empirical correction for the effect of the wake. The empirical 
factor K is given by Glauert as a function of (c/t), the 
values being derived by fitting equation (91) to the. experimental 
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data of reference 17. This m&e tern differs fundamentally from 
the wake correction of equation (67) in that the correction in 
this ca3e consists of a single tern which varies as (c/h); ■"■ 
whereas the correction in equation (67) comprises two terms, one 
of "which varies as (c/h) and one of which varies as (c/h) 2 . 
Equation (67) gives corrected results which agree as closely with 
the free-air values as do the results obtained with equation (91). 
It has the advantage that it is generally applicable to all air- 
foils and does not depend upon the experimental results of tests 
of specific sections. 

In suasrary ; the corrected data of figures 5 to 9 indicate, 
for the most part, that when the flow is ma.intai.ned strictly two- 
dimensional , the theoretical corrections for the tunnel- wall 
interference are, for low Mach numbers at least, accurate up to 
chord-height ratios of unity. The high accuracy observed at. the 
larger values of (c/h) must, however, be regarded as fortuitous 
since the theoretical analysis is predicated upon the assumption 
that the chord-height ratio is small enough that all points of 
the airfoil may be assumed to lie on the center line of the tunnel 
and that powers of (c/h) higher than the second may be neglected. 
It is thought that, at low Mach numbers, chord-height ratios as 
high as 0.7 are permissible if the tests are conducted only for 
the purpose of obtaining drag characteristics at low values of 
the lift coefficient. However, care must be exercised in 
ascertaining the maximum chord-height ratio permissible in any 
particular case to insure that the interference lift represented 
by equation (63) is not of such nature and magnitude as to affect 
the general character of the flow in the boundary layer along the 
surface of the model. In tests conducted to determine the aero- 
dynamic characteristics", of a model up to and beyond the maximum 
lift/ it is believed that the chord-height ratios must be kept 
to much lower values, At low Mach numbers, chord-height ratios 
up to 0.4 are probably permissible; however, there are no experi- 
mental data available at present to support this contention, 

• As noted previously, no experimental data could be found 
which would permit an evaluation of the accuracy of the calcu- 
lated effects of compressibility upon* the wall- interference 
corrections. Most certainly, as the test Mach numbers increase, 
the permissible chord-height ratios must decrease. Theory indi- 
cates that as long as the velocities induced at the position of 
the airfoil by the wall interference are small as compared with 
the velocity of the undisturbed stream, the corrections developed 
in this paper are valid even though the stream Mach number exceeds 
the critical for the airfoil under test. However, as previously 
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noted , at and for a range cf Kaon numbers below choking, the inter- 
ference velocities are no longer amali and the corrections are 
invalid.. The extent of this range is unknown. It should he 
emphasized that the flow pattern at and in the immediate vicinity 
of choking does not correspond to any flow pattern obtainable with 
the airfoil in free air; so the test results in this range cannot 
be corrected by any method. 

For zero Mach number- (i.e., for an incompressible fluid) , the 
results of the present paper can be compared with Goldstein's 
particular corrections for airfoils having small thickness and 
camber and small force coefficients. For an airfoil on the center 
line of the tunnel , equations (138), (139), (HO), (143), and 
(144) of reference 3, together with the expressions of appendix 
5 of reference 4, give the following equations for the velocity, 
angle of attack, and aerodynamic coefficients in an incompressible 



fluids 





(92) 



c d = c d ' \ 1 



i 



2 . 2 " 




, Here, the moment coefficients are for moments about the midchord, 
ana /c^A is the moment coefficient at zero lift in free air. 



The quantities C 0 and - C 2 * are determined by the shape of the 
base profile according to the equations 
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c« = 



c 2 



2 r h 

K J c sin G 



ae. 



/ cos 26 
/ c sin 6 



(93) 



In deriving these equations, the notation of references 3 and 4 
has been changed to agree with that of the present paper , and the 
lift-curve slope in free air assumed to have its .theoretical 
value of 2ic. 

The- corresponding corrections as obtained by setting 
M 1 - 0 in. equations (25), (62), (63), (64a), and (67) of the 
present paper are 



V = VM 



1 + Aa + Tcd'J • 



a = a' + {c z * + 4c nc « 



c z = oj! \ 1 - a - 2Act- 2Tc d ' 
L 

c d = c d' f 1 " 3/kT " 2Tc d' } 



c m-. - c m 



2Aa- 2Tc d 



J 



(94) 



The last of these equations is obtained from equations (.62) and 

1 

(63) by means of the relation = c m c ~ + -g C Z # 



The correction terms involving o in the two equations for 
the velocity are equivalent, except that the factor A, which . 
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appears in the equation of thy present paper ; is replaced in 
Goldstein's equation by the quantity (2C G - C 2 ) . Equation (93) 
gives 



- TT rr 

r9 o P n 4 / yt i - .cos ?,o Aa . .a / yt ■ . . 

o 2 KJ c sm6 ix f c 

o ' o 



-which becomes , after substitution from equation (44), 



i 



it ^ \ c / \ c/ c. 



(95) 



o - 



where ■ A is the cross- sectional area of the airfoil. The factor 
A can be expressed in analogous form by means of equation (13.05) 

1 2 

of reference 2. Since - Ac . is equivalent to the quantity 
in reference 2, this equation becomes 

fi A + A V / \ 

A =5 ^1 (96) 

5T C " 

where A v is the so-called" "virtual area" of the base profile. 

The virtual area of a given body in two-dimensional flow is de- 
fined as the area occupied by a fictitious quantity of fluid 
having a uniform density p and velocity V and possessing a 
kinetic energy equal to the total kinetic energy of the field of 
flow about the same body when it is moving forward with a steady 
velocity V through an unlimited expanse of incompressible fluid, 
of density p. The. magnitude of the virtual area depends upon 
the shape as well as upon the size of ^ the body, It is seen that 
the first correction term in the velocity equation of the present 
paper (which for the incompressible case is simply the result 
originally derived by Lock.) has a somewhat higher value than the ■ 
correction term of the Goldstein equation. The Goldstein equation 
contains no term corresponding to the term Tc^' in the equation 
of the present paper. Goldstein includes this correction, however, 
in the equation for the determination t of the stream velocity from 
measurements made at -the tunnel wall upstream, of the model. 



The Goldstein equation for the correction of the measured, 
drag coefficient. likewise differs from that of the present paper 
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"by the replacement of the factor A by the quantity . (2C 0 - C 2 ) 
and by the omission of the term in Tc^ 1 . If Goldstein's equation 

is applied to the experimental. results of Fage given in figures 
12 -and 13, it is found that there is. little to choose between the 
corrected results given by the two equations , except in the case 
of the circular cylinder where the results obtained from the equa- 
tion of the present paper are better. 

The corrections to the lift and moment coefficients as de- 
rived by Goldstein differ markedly from those of the present 
paper in that Goldstein's equations contain no terms corresponding 
to the 2Ag and 2^o^ terms which appear in the equations of 

this paper. As has been noted previo ; asly, the . 2Tc ( j t term. is 
accounted for indirectly in the determination of the apparent 
stream velocity. A term ef the type 2Ao is necessary, however, 
to correct the measured coefficients for the increase in dynamic 
pressure caused by the interference between the walls and the air- 
foil thickness. 

, Since the moment coefficient at zero lift is the same about 
any axis and since the change from the free air to the measured 
moment coefficient in the correction to the angle of attack will 
introduce only differences of the second order in a, Goldstein's 
equation for the corrected angle of attack may be written with 
sufficient accuracy as 



In this equation, the part of the correction due to the moment on 
the airfoil is .constant, its value depending only upon c^' for 



zero. lift; whereas in the corresponding equation of the present 
paper the part of the correction due to the moment varies with 
the angle of attack. . This difference • is of small consequence in 
most applications; however, the equation of the present paper, 
which includes, the actual variation in moment, may be somewhat 
the more accurate, especially at high angles of attack. 

The compressibility factors which appear in the complete 
equations of the present paper are comparable with the results of 
Goldstein and Young (reference 5). The equation for drag as given 
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in reference 5, when expressed in coefficient form and altered to 
agree with the notation of the present paper,. can he written 



c d = c d » J 



[l-(M') 2 ] : 



Ac 



[l-fM') 3 ]. 



(97) 



This equation is obtained by modifying equation (91) to include 
the effect of compressibility. Comparison of the compressibility 
modifications of equation (97) with those of the corresponding 
terms of equation (67) reveals that the compressibility factors 
appearing in the first correction terms differ by the inclusion 
of a term -(M r ) s in the numerator in equation (67). This 
difference arises from a failure to note in the development of 
equation (97) that in a compressible fluid the dynamic pressure 
in the tunnel is affected by the change in density which ac- 
companies the change in axial velocity. The compressibility 
factor of the 3eccnd (or wake- correction) term of equation (97) 
is not comparable with the compressibility factors of the wake- 
correction terms of equation (67) because of the fundamental 
difference in the nature of the corrections already pointed out 
in the discussion of equation (91). The compressibility factors 
in the equations for lift, moment, and angle of attack in reference 
5 agree with those appearing in the corresponding terms of the 
equations of the present paper. It should be noted, however, 
that the lift and moment equations of reference 5 include no 
corrections for the difference between the true and apparent 
dynamic pressures in the tunnel. 



CONCLUSIONS 

Airfoil data obtained from tests in a two-dimensional-flow 
wind tunnel can be corrected to free-air conditions by means of 
the following equations: 



1 + 



[l.(M0T° /2 1 - (MO 2 
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(63) 



(67) 



where 

hod 

* 2 /cV 

and A is a dimensionless factor the value of which depends upon 
the shape of the "base profile of the airfoil. (See equation (3) 
and tahle I.) The remaining symbols are defined in appendix B. 
Numerical values of the functions of M 1 which appear in these 
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equations are given in table II, Experimental pressure distri- 
butions also can be corrected by a method "Outlined in the text. 

The corrections derived should be valid theoretically up to 
a Mach number near the choking value, which is the maximum Maoh 
number attainable in the wind tunnel. The choking Mach number is 
shown to be the stream Mach number at which a Mach number of unity 
is attained locally across the tunnel either (1) at the position 
of the airfoil because of the reduction of the available flow 
area occasioned by the presence- of the airfoil, or (2) downstream 
of the airfoil as a result of the influence of -the . airfoil drag , 
upon the flow in the wake. The choking Mach number can be esti- 
mated by means of equations presented in the report. 

Insofar as can be ascertained from the small amount of 
experimental data available, the correction equations are appli- 
cable at low Mach numbers for values of the chord-height ratio 
(c/h) as high as 0,7 if the tests are conducted for the purpose 
of obtaining drag characteristics at low values of the lift co- 
efficient. In tests conducted to determine the aerodynamic 
characteristics of an airfoil up to and beyond the maximum lift, 
it is thought that a chord-height ratio of 0.4 is permissible at' 
low Mach numbers, although there is no experimental evidence 'to 
support 'this contention at present. At high Mach numbers the 
permissible chord-height ratios must logically be expected to. 
decrease. In particular, if the critical speed is exceeded, it 
is probable that only very small values, of (c/h) are permissible. 
There are at this time no experimental .data available on this . 
aspect of the problem. 

Comparison of the results of' the present paper with those 
of references 3, 4, and 5 reveals certain differences as noted 
in the section Discussion. 
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APPENDIX A 



THE VELOCITY FIELD FOE A SOURCE AND FOR A SYMMETRICAL AIRFOIL 

IN A COMPRESSIBLE STREAM 

The velocity induced at a point in a compressible fluid 
stream by a single fluid source can be found to a first degree of 
approTl.Tna.ttoh by a modification of the method used by Glauert 
(reference 7) for the consideration of a vortex in a compres- 
sible fluid. To this end, a system of polar coordinates is intro- 
duced. The origin is located at the source and the polar axis 
extends; downstream parallel to the velocity V of the undisturbed, 
stream. (See fig. 14.) The resultant velocity U at any point 
L(r, c P) is defined by the velocity components w and n parallel 
and normal, respectively, to the radius vector. 

The condition for irrotational motion requires that at all 
points in the field 

dCrn) . = o (Al) 

3r an> 

The equation of continuity is 

W™) - + ^ + ™5£l+IL^ = o (az) 
Br &p Pl br Pz aq> 

where pj is the density of the fluid at any point. 

The source strength (mass flow per unit time) is denoted by 
Q. Then, for any circle enclosing the source, considerations of 
symmetry and continuity, respectively, provide the two integral 
relationships 

2TT 

^ nr dS> = 0 (A3) 
and 0 2tt 

p z wr d9 = Q (A4) 

o 
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The radial and circumferential components of the velocity 
may "be expanded in the series 



w = V j^cos 



CO V 

3=1 / 



n = -V 



7 



where A g and B s are functions of <P. If r is large, it is 

sufficient to retain only the first terms of each power series - 
so that ; 



w = v/cos<P + £^- 



(A5) 



J L . 

n = -V ^sin <+> + | j J 

To' the first power in (l/r), the square of the resultant 
velocity is 

(?) = ^~ 2 — = 1 + f (A cos <P + B sin CP) 

For reversible adiabatic flow, the local density p x is related 
to the density p of the undisturbed stream by 



! 7- l 



J 



where M is the Mach number of the undisturbed stream and y is 
the ratio of the specific heats. Thus, to the first power in 
(l/r), 



1 - — (A cos <¥> + B sin <f) , 



(A6) 
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The solution is now obtained "by inserting values from 
expressions (A5) and (A6) into the fundamental equations. Sub- 
stitution of (A5) into the equation for irrotational motion (Al) 
requires that A shall he a constant. Substitution of (A5) and 
(A6) into the equation of continuity then gives 

|| (1 - M 2 sin 2 cp) = M 2 (A cos 2<P+ B sin'2<P) 

which becomes upon integration 

B(l - M 2 sin 2 <P) = - M 2 A sin 29 + C (A7) 

2 

where C is a constant. The integral equations (A3) and (A4) 
become, respectively, 

BdCp« 0 (A8) 



and 

2TT 



B sin 2<P d<P = 2rtA (1 - - JL ' 4 ■ ' (A9) 

V 2 J PV 



o 



Substitution of the expression for B from equation (A7) 
into the integral equation (A9) gives 



A = -a- 
2npV 



while substitution into equation (A8) showB that C => 0. Thus, 
from equation (A7), 



4itpV 



f M 2 sin 2<? 1 
(1 - M 2 sin 2 9)Vl - M J 



The expressions for the velocity components therefore become 
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w = V cos <P + 



&tpr - M 2 



n = -V sin <r> - _S_ — 

2itpr 



M 2 sin<P cos V 



M 2 sin 2 cp) /y/T"- 



(A10) 



For a Mach number of zero these equations reduce to the vell- 
Icnovn results for a source in an incompressible fluid. 

From equations (A10) the velocity components . u and v, 
parallel and perpendicular, respectively, to the direction of the 
undisturbed stream, are found to be 



u = V + 



v = 



2rtpr 

_iL 

2rtpr 



cos 



CO 



*Jl - M 2 (i - M 2 sin 2 r P) 



\ fl - M 2 sin 9 
~" 1 - M 2 sin 2< P 



(All) 



The drag experienced by the source can "be determined by- 
evaluating the integral 

2TT 




Vl cos C P+ pjw(v cos ^ - n sin ^ > rd <P 



D = 



over any circle enclosing the source. To the accuracy previously 
employed, the pressure at any point is 



p, = p - 2^ (A cos <P + B sin cp) 
I r 



(A12) 



Insertion of this expression, together with (A5) and (A6)) into 
the equation for drag gives finally 



D = - VQ 



(A13) 



which is the same as for a source in an incompressible fluid. It 
is apparent from considerations of symmetry that the lift force 
of the source is zero'. s ^ 
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The results of equations (All) can be used to study the field 
of flow about a symmetrical airfoil at zero angle of attack in a 
uniform stream. Such an airfoil can be represented by a suitable 
system of sources and sinks distributed continuously along. the 
chord line. If the notation of figure 15 is used, the vertical 
velocity v^ induced in an incompressible stream at a given point 

( x oj yo) on the surface of the airfoil is 

T1 . _L f*** ("SU . J, f_^_ a («*\ a, (au) 

acp J r \dx / atp J (x 0 -x) fc +7o 8 \te I 

o o 
* /dQA 

where ( 1 j is the strength per unit length of the source-sink 

distribution in an incompressible stream. From the second of 
equations (All) it follows that the velocity v c ' at the same point 
in a compressible stream is 

-JL / P ^/^T^gincp ^^ dx 

dx / 




r(l-M 2 sin 2 qp) Vdx 



or 



!i=5f f y o „ ( fbft dx (A15 ) 

*P J (x 0 -x) 2 + (1-M )y 0 S vdx / 



\/l-M a 
2*0 ^ 

o 

where is strength of the source-sink distribution in 

the compressible case. For any given airfoil of small thickness 
the condition that the flow shall be tangential to the surface 
of the airfoil requires that v c = v± at all points on the 

surface. This fact can be used to relate the source- sink distri- 
butions for a thin airfoil in the compressible and incompressible 
streams by considering the limiting forms of equations (A14) and 
(A15) as y Q approaches zero. 

Consider first the limiting form of .equation (A14), which 
may be written 

. ' . T1 .U» *2- f 1 (iSft te (AM, 
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It is. seen; that even for y 0 = 0. the integral in this equation 

is finite when evaluated ovor any interval of integration not 
including the. point .x = x 0 . In the limit, the contribution of 

■such intervals to the right-hand side of the equation is there- 
fore zero, and the equation may be written ' 



lim , lim 
€ -r^O yo->0 




(x 0 - x) 2 + y 0 2 \&x J 



dx 



In evaluating the limit in this equation, care must be taken that 
the limit with respect to y 0 is taken first in every case. 
Integration "by parts gives 



v i " 



lim 



lim 



e— >0 ■ y 0 — >0 



(tan -1 ) 




y 0 



x 0 +€ 



x o" € 



(tan" 1 ) dx 

y 0 



By virtue of the first mean value theorem for integrals 



(A17) 



reference 



18, p. 65) the integral term in this equation may be written 



x 0 +e- 



X, 



.x^-x. 



x o" € 



(tan- 1 ) 0 dx=. 

y 0 . 



'd 2 Q 
V dx 



_2 /x= 



r 



(tanf 1 )* 0 "* dx 

y 0 



x 0 -e 




-1 X 0" X 

(tan ) dx 



where (x 0 -e) % i x = x 0 and x 0 = £ 2 ^ (x 0 +c). The division into 
two integrals is necessary to ensure that the conditions under 
which the mean-value theorem is applicable are fulfilled; namely, 
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that 



PS 



is continuous and that (tan" 1 ) — ^ 'has the 



same sign throughout the interval of integration. Integration 
gives 



z 0 + € 



dx = 



Le" (*x 2 ) x= i 



X 



e (tan ) £ log 

2. 



(■•£■)] 



In the limit , the valuS. of the terms in the second, "bracket in 
zero. Thus, only the first term need be retained in equation 
(A17), which' may now be written 



= lim 



lim -Lftan" 1 ) -1 



/dQi\ 
J x=: 



+ (~) 
X=X Q +€ \** J X=X Q -€ 



lim 



iim (tan" 1 ) 
— > n y o 



Thus the limiting form of equation (A14) "becomes finally 



v,. and I ) 

1 - WJ 



2p \dx J 



(A18) 



where v,. and may now "be considered as certaining to the 

same general- chord wise station x. 
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The limiting form of equation (A15) can "be. found in similar 
fashion. In this case integration "by parts gives in place 'of 
equation (A17) 



v c = lim 



lim 

o y^o 



2itp 



\dx 



—(tan" 1 )— 2_ : 

M 2 V 0 -%P^' C 



x o +£ 



x o" s 



(tan' 1 )- 0 



dx 



(A19) 



As before, the value of the integral term in this equation is 
zero. The limiting form of equation (A15) becomes finally 



2p \dx j 



(A20) 



which is the same as (A18). 



Since for any given airfoil v c = Vj^ at all chordwise 
stations, it follows from (A18) and (A20) that 



(A2l) 



that is, the source-sink distributions necessary to represent 
any given thin symmetrical airfoil in a uniform stream are 
identical for the compressible and incompressible case. 

This result can be used to calculate the effect of com- 
pressibility upon the field of induced .velocities at a consider- 
able distance from the airfoil. The increase in longitudinal 
velocity at a large distance y 1 directly above or below the 

midchord point of the airfoil' in an incompressible fluid is 
approximately 



dz 
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By virtue of the first of equations (All), the corresponding 
velocity at the same point in a compressible fluid is approximately 

c ' 

(u c - V) - L_ f (2 - x\ P*) dx 

o 

Thu3, in view of equation (A21), 

(u. - V) = 1 T (Ui - V) (A22) 

that is, in a compressible fluid the increase in longitudinal 
velocity at a point a considerable distance directly above or 

below a symmetrical airfoil is l/(l - M 2 ) 3/a . times the increase 
in longitudinal velocity at the same point in an incompressible 
fluid. 

The foregoing results can be used, also to determine the 
effect of compressibility upon the drag of an airfoil in a stream 
having a longitudinal pressure gradient. Consider an undisturbed 
nonuniform stream having at some given point a velocity V, a 
density p, and a streamwise pressure gradient dp/dx. By 
virtue of Bernoulli's equation, there must be at this point a 
velocity gradient' 



4S « -L & , (A23) 

dx pV dx 

This holds true both in the compressible and the incompressible 
case. The velocity v a small distance x from the point in 
question is then 

v ,V + x^ = V- ^^ (A24) 
dx pV dx 

As a result of equation (A13), the drag experienced by an airfoil 
placed at this point in the stream is, for both the compressible 
and incompressible cases, 
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where (dQ/dx) is, as before , the strength per unit length of 
the source-Gink system necessary to represent the airfoil. In 
order to fulfill the condition that the airfoil is a closed 
body., the source- sink system must he such that 



c 

r 



o 

Thus the drag is fir-ally 



/ (!) 



dx = 0 



pV dx, J 



m 



D = 4 : $2 / X dx (A25) 



in both tiie compressible and the incompressible cases. If the 
streamwise pressure gradient is small, equation (A21) is still 
applicable; that is, the' source-sink distributions necessary to 
represent the airfoil in the compressible and incompressible 
cases are identical. It therefore follows from equation (A25) 
that the drag of an airfoil in a, stream having a longitudinal 
pressure gradient is unaffected by fluid compressibility. 



\ 
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APPENDIX B 
LIST OF IMPORTANT SYMBOLS 

c airfoil chord 

t airfoil thickness 

h tunnel height 

A a factor depending upon shape of base profile (see equation 
(3) and table I) 

G* ?^ factor depending upon size of airfoil relative 

• to tunnel 

T - (2\ ; factor depending upon size of airfoil relative 
to tunnel 

a angle of attack ■ 

Cj section lift coefficient 

c mc section quarter- chord-moment coefficient 

c m c section m\& chord-moment coefficient 

2 

c,j section drag coefficient 

V stream velocity 

M Mach number 

M'c^ apparent Mach number at choking 

R Reynolds number 
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7 — ; ratio of specific heat of gas at constant pressure to 

to specific heat at constant volume (for air y = 1.4) 
H total head 
p static pressure 
q dynamic pressure 
p mass density 
ji coefficient of viscosity 
T absolute temperature 
Y c speed of sound 

1+T] . compressibility factor (see equation (71) and fig. 4) 

D . .section drag 

AD section drag due to stresunwise pressure gradient 

P chordwise lift distribution in coefficient form 

P e interference lift distribution (see equation (78) and 
table III) 

P2 local pressure coefficient (see equation (69)) 

local pressure coefficient (see equation (68)} 

x coordinate of points on chord line as measured from 
leading edge 

6 angular coordinate of points on chord line (see equation 
(43)) 

r radial distance in polar coordinates 

cp polar angle in polar coordinates (positive counter- 
clockwise) 

t^ projected thickness of airfoil 
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t e effective thickness of airfoil 



y t ordinate of base profile 
slope of mean- camber line 



dx 



Q source strength 

vorticity per unit length of chord line 



dx 



u horizontal component of velocity 
v vertical component. of velocity 

n circumferential component of velocity in polar coordi- 
nates (positive counterclockwise) 

w radial component of velocity in polar coordinates 

A geometrical area of airfoil section 

A v virtual area of airfoil section 

A ? cross- sectional area of empty tunnel 

A_ • minimum cross- sectional area between model and tunnel walls 
m 

"Aj local cross- sectional area of stream tube 
A n Fourier coefficients (see equations (45) and. (49)) 
Superscripts. 

( J ) when pertaining to fluid properties, denotes values exist- 
ing in' tunnel far upstream from model; when pertaining 
to airfoil characteristics, denotes values in tunnel, 
coefficients being referred to apparent dynamic pressure 
q« ' 

(") denotes fluid properties far downstream from model 

(*) denotes airfoil characteristics in tunnel as coefficients 
referred to true dynamic pressure q 
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Subscripts 

c denotes values in compressible fluid (excepting V c ) 

i denotes values in incompressible fluid 

I denotes local conditions at point in fluid 

s denotes conditions existing far downstream vhen airfoil 
and wake are replaced by source 

m denotes conditions at minimum, cross- sectional area between 
airfoil and tunnel walls 

L - denotes values on lower surface of airfoil 

U denotes values on upper surface of airfoil. 
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TABLE III 

VALUES OF ? e AT S^WQAHD CEOHD'/ixSE STATIONS 



1 

z/c 




x/c 




0 


0 


0,40 


1,2475 


0,005 


0.17S6 j 


0.45 


1.2669 


0.0075 


0.2198 


0.50 


1.2732 


C.0125 


0.2830 


0.55 


1.2569 


0.025 


0,3976 


0.60 


1.2475 


0.050 


0.5550 


0.65 


1.2146 


.0.075 


0.67,07 


0.70 


1.1670 


0.10 


0.7639 


0.75 


1.1027 - 


0.15 


0.9093 


0.80 


1.0186 


0.20 


1.0186 


0.85 


0.9093 


0.25 


1.1027 


0,90 


0.7639 


0.30 


1.1670 


0.95 


• 0.5550 


0,35 


1.2146 • 


1.00 


0 
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Figure 1.- Source system for analysis of wake effect. 
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Figure 2.- Meaii-camfcer line in free air. 
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Figure 3.- Mean-camber line in tunnel. 
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Mach number, M 



Figure 4.~ Compressibility factor. 
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Figure 6.- Lines of sonic speed at the position of the airfoil 
after choking. 
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(a) Uncorrected for tunnel-wall interference 



. Figures 9a, b.- Lift characteristics for KACA 0012 airfoil section. 
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(b) Corrected for tunnel-wall interference 

figure 9b. 
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(a) Uncorrected for tunnel-wall interference 

Figures lOa.b.- Lift characteristics for NACA 23012 airfoil section 
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(b) Corrected for tunnel-wall interference 

Figure 10b. 
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(a) Uncorrected for tunnel-wall interference 

Figures 11a, b.- Moment characteristics for NACA 0012 airfoil 
section. 
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(b) Corrected for tunnel-wall interference 

Figure lib. 
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0 Corrected for wall interference, c 
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(Equation 67) 
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Figure 12.- Profile drag for three symmetrical Joukowski air- 
foils at zero angle of attack. 
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Figure 13*- Profile drag for three symmetrical 
bodies at zero angle of attack* 



